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Abstract 

It is classical that well-known identities and properties of partial quotients fur¬ 
nish rational approximation in F[p —1 ]]. For a rational function, this is the extended 
Euclidean algorithm in F[x], Berlekamp’s heuristic solution of the ’key equation’ 
essentially approximates an element of F[x] with constant term 1 via a quotient of 
reciprocals, and his solutions satisfy a number of identities. In earlier papers we 
gave a solution of an analogous problem using D[x -1 ,x], D a commutative domain. 

The linear complexity (of a finite initial subsequence) of an infinite sequence 
over F has been related to the degrees of its partial quotients by Mills, Cheng, 
Niederreiter and others. We use first principles and induction to relate these linear 
complexities to the degrees of its partial quotients. 

Berlekamp has also described the set of solutions of the key equation. We de¬ 
fine a pairing of minimal solutions and a ’minimal system’ of a finite sequence 
over D. Examples are classical approximation in F[[x -1 ]] and approximation using 
D[x -1 ,x]. We use minimal systems to generalise results of Massey and Niederre¬ 
iter to arbitrary solutions, including numerators. This includes explicit and unique 
decomposition of both parts of a solution into a sum of (polynomial) multiples of 
solutions with minimal degree denominators. The unique multipliers also satisfy 
degree constraints. 

We give several applications to gcd’s of sequence polynomials and relate partial- 
quotient solutions to solutions derived using F[x -1 ,x]. We give a precise count 
of the number of solutions when the field is finite. Our final application concerns 
when the first component of a minimal solution vanishes at some scalar; a simple 
modification of our approach gives a new solution, the first component of which 
does not vanish at the scalar and which has minimal degree. We also describe the 
corresponding set of solutions. This simplifies and generalises work of Salagean. 

We conclude that numerators (or second components) of solutions can play a 
significant role in proofs of properties of denominators (or first components) and 
that they enjoy similar properties. 

Keywords Berlekamp-Massey algorithm, continued fraction, key equation, Laurent 
series, linear recurrence, minimal polynomial, partial quotient, rational function. 

*School of Mathematics and Physics, University of Queensland, Brisbane, Queensland 4072, Australia 
(email: ghn@maths.uq.edu.au). 


1 



Contents 


1 Introduction 

1.1 Background. 

1.2 Overview. 

1.3 In More Detail . 

1.4 Standard Notation. 

1.5 Guide to Additional Notation 


3 

3 

3 

4 
6 
7 


2 Sequence Basics 

2.1 Rational Approximation and Solutions . 

2.2 Linear Recurring Sequences. 

2.3 Annihilating Polynomials and Solutions . 

2.4 First Examples and the Key Lemma [2.61 

2.5 Linear Complexity and Minimal Solutions 


8 

8 

8 

_2 

10 

li 


3 Minimal Solutions via Partial Quotients 

3.1 Continued Fractions . 

3.2 Geometric Sequences. 

3.3 Essential Sequences and the General Case 


0 

12 

15 

15 


4 An Inductive Construction of Minimal Solutions 

4.1 Pseudo-Geometric Sequences. 

4.2 Essential Sequences or, a Tale of Two Triples. 

4.3 The Inductive Theorem and the Corresponding Algorithm 

4.4 A Worst-Case Analysis. 

4.5 An Identity for /i and /i'. 


18 

18 

19 

20 

24 

25 


5 Decomposition 

5.1 A Pairing. 

5.2 Geometric Sequences. II .... 

5.2.1 Annihilating Polynomials 

5.2.2 Solutions. 

5.3 Essential Sequences. II. 

5.3.1 Annihilating Polynomials 

5.3.2 Solutions. 



6 Some Applications of Decomposition 

6.1 Sequences over a Field. 

6.2 Non-Vanishing Annihilating Polynomials 


[32 

32 

33 


2 




















































1 Introduction 


1.1 Background 

Let F be a field. Approximating the generating function Six -1 + S^x -2 + • • • G x -1 F[[x -1 ]] 
by rational functions ($/$ (where i > 0) is well-known; q[ l \ are known as its partial 
quotients or rational convergents. An important identity is 

- 9!‘ ) 9?- 1) = (-D- 1 . (1) 

Obtaining partial quotients uses division in the field of Laurent series in x _1 , written 
F((x -1 )). When the above sum is a rational function, this is the extended Euclidean 
algorithm. See also [8] for connections with linear recurring sequences. 

A second example is Berlekamp’s iterative solutions cuW/cr^ of the ’key equation’, 
where 0 < i < n, [U Section 7]. (The integer n is related to a decoding problem.) It 
is essentially rational approximation of 1 + Six + • • • + s n x n € F[x] using reciprocals of 
polynomials. It uses ’auxiliary solutions’ yh)/ r M which satisfy = x®, [TJ 

Theorem 7.42], The set of solutions was discussed in [T, Theorems 7.43, 7.44], 

A simplification of Berlekamp’s algorithm appeared in [7} Algorithm 1], This interprets 
aO) as a ’connection polynomial of a minimal-length linear-feedback shift register (LFSR) 
which generates s — Si ,..., s n \ It is known as the Berlekamp-Massey algorithm. The 
minimal length is called the ’linear complexity’ L n of s. The set of connection polynomials 
for all LFSR’s of length L n which generate s was given in pj, Theorem 3]. 

Connections between these two types of rational approximation e.g. between the linear 
complexity of Si,..., S n and the degrees of the denominators q'f 1 have been discussed in 
0 . i and [15], which depend on [lj. In [9] Theorem 1] this was done independently of 
[lj and [7]. 

A third example appeared in m Our goal was a faithful redevelopment and exten¬ 
sion of [I]; we were unaware of [9] and Macaulay’s inverse systems, see e.g. [TO] . We 
discussed rational approximation of s 0 + • • • + Si_ n x 1-n G D[x _1 ] using Laurent polyno¬ 
mials D[x~ x , x], where D is a commutative domain. We write our solution as fi — (// 1 , /i 2 ) 
and call /v. 1 a ’minimal polynomial’ of s. The linear complexity L n of s is the degree of /j,i 
and the reciprocal of fii is a connection polynomial of an LFSR generating s, [IT]. When 
D is a field, our approach has applications to the above decoding problem and to control 
theory, see for instance [12, Section 8] and pT] Example 4.9]. 

1.2 Overview 

Our overall goal is to unify and extend some results in the literature related to the rational 
approximation of generating functions of infinite and finite sequences. In our approach, 
numerators play a significant role. 

We revisit [9] Theorem 1], which has two parts. We give an inductive proof of the first 
part on linear complexity and partial quotients. Our proof is from first principles, using 
the basic definitions for finite sequences from m We also prove the converse. 
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We also derive an analogue of Identity (00) for our minimal solutions, Proposition 14. 201 
This enables us to ’decompose’ solutions and determine the set of all solutions for a finite 
sequence over D whenever we have a ’minimal system’ for s. Partial quotients (with 
D = F) also provide a minimal system. In this way we generalise the second part of [9( 
Theorem 1] to all solutions. We conclude with some applications of decomposition. 

Note to the reader: we consider the partial quotients for an infinite sequence over a 
field only; we have not extended [2J and [9] to commutative domains. In some situations, 
we apply our results to finite sequences over a field e.g. Proposition 12.81 Proposition 14.141 
on monic minimal polynomials and Corollaries 16.11 - 16.41 We have included a number of 
examples; some reappear intentionally in different guises as an expository aid and others 
are inductive bases for later theorems. 

1.3 In More Detail 

We begin with basic concepts for infinite sequences over F, denoted So, S- 1 ,... and finite 
sequences over D, denoted so,..., si_ n where n > 1; this indexing agrees with Macaulay’s 
inverse systems in pD] and with finite sequences in m 

We can regard finite sequences as trivial (s = 0,..., 0), geometric or ’essential’, Propo¬ 
sition [2JJU Geometric sequences are those of high school, defined by s 0 ^ 0 and a common 
ratio. Equivalently, they satisfy L n = ■ • ■ = 'Ll = 1. ’Essential’ sequences on the other 
hand satisfy L n > L 1 > 0 and predominate: geometric sequences may become essential 
on adding a term, but never the reverse, see Proposition 12.111 We summarise this using 
a state diagram (T is the start state, ’G’ denotes ’geometric’ and ’E’ denotes ’essential’; 
we have suppressed transitions between the same state): 



n> 2, Ai^O 


where s ^ 0,...,0 and denotes a ’discrepancy’. Unfortunately this subdivision of 
sequences does not appear in [TJ , which renders the Berlekamp-Massey algorithm harder 
to understand. If n > 2 and s — sq, ... , s x _ n is essential then 

n' = max{l < j < n : Lj < L n } 

is a well-defined integer, 1 < n'< n, and we have the important subsequence s o, • • •, si_ n '- 
We next discuss partial quotients (when S ^ 0, 0,...) as these are classical and less 
detailed, being based on division in F((x -1 )). First we treat the base cases in Propositions 
13.4113.61 We obtain an inductive proof of the first part of [9J Theorem 1] and its converse, 
Theorem 13.71 This gives a similar state diagram for S, Proposition 13.81 
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Then we revisit HU, restricting to geometric and essential sequences over D only. This 
new approach is simpler, see Theorem 14.91 the corresponding Algorithm 14.121 is valid for 
all sequences and is virtually identical to HD Algorithm 4.6]; we compute p = (pi,p 2 ) 
rather than (pi,a:p 2 ). (Apart from Lemmas 14.21 and 14.71 — which the interested reader 
may verify — this paper is independent of HUO 

For sequences over a field, there is a ’normalised Algorithm 14.121 which computes a 
monic p 1; Proposition 14.141 This has been implemented in COCOA, [3]. We also prove 
an identity for the sum of the linear complexities of s. This seems to be new and gives a 
simple analysis of Algorithm 14.121 see Proposition 14.171 

We use Algorithm 14. 121 to define an element V s £ D \ {0} and prove the identity 

H-2 n[ - pi p 2 = V s (2) 

where p 7 is either (1, 0) or a minimal solution for s o, • • ■, Si_ n '. This is our analogue of 
Identity ([!]) and the identity of Berlekamp mentioned above. Identity ([2]) easily implies 
that for any / = (/i, / 2 ) E D[x] 2 

V. 5 ./i = rn // 1 - rn p] (3) 

where m = / 2 pi — fi p 2 and m! = / 2 p] — /ip 2 . This is a special case of a pairing 
D[x] 2 x D[a;] 2 —> D[rc] defined by p and p'. 

In fact, essential sequences exhibit a ’minimal system’, a stronger property than (|3l) . 
Definition 15.81 We show that if / is a solution and we have a minimal system, then these 
multipliers (i) satisfy degree bounds and (ii) are unique when the degree of fi is at most 
n; in this case we call Q a ’decomposition’ of V s f\. And V / 2 satisfies a similar identity 
with the same multipliers i.e. we have a decomposition of V s /. This yields the required 
description of all solutions when we have a minimal system, see Corollary 15.151 Partial 
quotients also exhibit a minimal system. In this way, we generalise the second part of 0 
Theorem 1], It also strengthens [TT| Theorem 4.17] and has a simpler proof. 

There are lacunae for geometric sequences as they do not have a minimal system. 
However this does not embarrass us, as using ((/[°\ q^) = (p(, p 2 ) = (1,0) enables us to 
give alternative proofs in both the partial quotient and finite sequence contexts. Secondly, 
over a domain D, we have to work with ’pseudo-geometric’ sequences as the leading coef¬ 
ficient of p! may not be a unit of D. As these sequences are inherently simpler and easier 
to treat than the essential ones, we always discuss them first. 

We have included some applications. We show that for a sequence s over a field and 
solution (/i,/ 2 ) such that the degree of fi is at most n — L n (i) gcd(/i,/ 2 ) = 1 implies 
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that / is a minimal solution and (ii) the multipliers of Identity (J3J) satisfy gcd(m, m') = 
gcd(/i, / 2 ). We apply (i) to linear recurring sequences. We relate our minimal polynomials 
and partial quotients, Corollary 16.21 We also give a precise count of the number of 
solutions when |F| < oo. 

For our final application, we revisit some work of Salagean, H3J. Let a G D be 
arbitrary and suppose that n\ (a) = 0. We show that Identity ([3]) implies that the lower 
bound for the degree of an annihilating polynomial of s which does not vanish at a is 
M = L n + max{n + 1 — 2L n , 0}. We exhibit a solution of minimal degree 


- the polynomial x M /q — fi\ does not vanish at a by Identity ([3]). Algorithm 16.81 is a 
one-line extension of Algorithm 14.121 and is simpler than [lfU Algorithm 3.2], We also 
derive the corresponding numerator. In fact the bound in Theorem 16.71 and the set of 
minimal polynomials in Corollary 16.121 were stated without proof in [TT] and used to jus¬ 
tify Algorithm 3.2, loc. cit. 

We thank the anonymous referee for a number of useful comments and suggestions 
which improved the presentation, and also the members of Projet Secret at INRIA, Roc- 
quencourt for their hospitality. 

1.4 Standard Notation 

For any set E containing 0, E* = E \ {0} so that N x ={1,2...}. As usual, = 0. 

Throughout the paper, D is a commutative domain with 1^0 and R = D[x]. For 
any a G D x and A C D, aA = {ab:bE A}. For / G R, |/| is the degree of / G R, with 
101 = — oo; the usual rules for arithmetic involving — oo apply. If / G R x , lc(/) is the 
leading coefficient of /. We often write / = x k g + h if /(x) = x k g(x) + h(x) where fcGN 
and g,h G R. For /, g G R, their product is written / g and we regard R 2 as an R-module 
via f(g,h) = (fgjh). 

A non-zero formal negative Laurent series over D is L = 'ff^ 00<i< i : LiX l G D((a;^ 1 )) 
where k G Z, Lj G D and L/, 0; we write 

k 

v(L) = k and [L\ = Li X 1 G x R 

i —1 

i.e. v : D((x -1 )) —> {—oo} U Z is the exponential valuation , with v(L) = — oo L = 0; 
v coincides with | | on R. It is elementary that v(L • M ) = v(L) + v(M), v(L + M ) < 
max{v(L),v(M)} and v(L + M) = max{v(L), v(Af)} if v(L) v(M). We also write v 
for the restriction of v to Dfa; -1 ] C D((x -1 )). We regard D((a: _1 )) as D[[a: _1 ]] + xR and 
use • for multiplication in D((x“ 1 )). 

We denote an arbitrary held by F. For continued fractions in F^x^ 1 ]] we use Pol(L) = 
Y^i=o LiX 1 G F[x]. As usual, F(x) C F((x^ 1 )) is the subheld of rational functions over F. 
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1.5 Guide to Additional Notation 

We include a table of additional symbols used in the paper to aid the reader. 


Symbol 

Meaning 

Q n 

sequence of n zeroes 

a, b, c 

elements of D 

cii 

[di, a 2 ,...] is the continued fraction expansion of S and z > 1 

Ann(s) 

set of annihilating polynomials of s 

bi 

Laurent series in z th iteration of partial quotient algorithm, i > 0 

6 5 6 S 

n + 1 — 2L n 

f,9 

elements (/i,/ 2 ), ( 91 , 5 - 2 ) of R 2 

h 

[A ■ s]/® 

(f,g) 

,/ 2 ,9i - /192 e R, the pairing of / and g 

ld s 

ideal of characteristic polynomials of S 

L, L(s), L n 

linear complexity of s — s 0 , ■ ■ ■, s i„ 

L',L( S '),W 

linear complexity of s', s essential 

m, m' 

(/,b), (/,b') respectively 

MP(s) 

set of minimal polynomials of s 

n 

a strictly positive integer 

n' 

the strictly positive integer maxi<j <n {j : Lj < L n }, s essential 

rii 

strictly positive integer \q[' 1-1 1 + g}*" 1 or 00 

V 

V, 5 (non-zero product of discrepancies) or (—1)* 1 

qd) 

z th partial quotient (q^, q^), z > 0 

s 

finite sequence so,..., si n over D 

n 

So 5 • • • 7 ^l— S—n 

s 

So + s iX 1 H-h s 1 n x 1 ~ n 

s' 

the subsequence s 0 ,..., s 1 „/ of essential s 

S 

inhnite sequence over F 

S 

Ao + S- yx- 1 + • • • 

S\n 

So, . . . , Si- n 

T s 

triple (n,//i, Ai) for s 

K 

triple {n', , A)) for s', s essential 

Ai = A(/ii; s,s n ) 

next discrepancy of // 1 e MP(s) 

a; 

next discrepancy of //( 

A 

either /i or as used in Section [5] 

b 

minimal solution (/Zu/Zo'l for s from Theorem |4.9| 

bi 

minimal polvnomial of s from Theorem |4.9| 

bi 

1 if s is pseudo-geometric, or minimal polynomial of s' 

V 

new minimal solution obtained from fi, fi' 

t 

solution with 0 (a) 0 0 constructed in Section 6, a G D 

(f), (j)', 0 

elements of R. 
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2 Sequence Basics 

2.1 Rational Approximation and Solutions 

Given an infinite sequence S = So, S_ i,... over F, rational approximation of the gener¬ 
ating function of S and continued fractions is classical. Consider the following problem: 
for n > 1, find a rational function x f 2 /fi G ¥(x) with |/ 2 | < |/i| such that 


(x f 2 /fi)i = Si for 1 - n < i < 0 (4) 

and |/i| is minimal. Let S_ = J2i<o^iX l G F[[x -1 ]] be the generating function of S. We 
can rephrase (J4j) as: find xf 2 /fi such that v(S_ — xf 2 /fi) < — n and x f 2 — [fi ■ S). 
Multiplying by fi, we equivalently require (fi, / 2 ) such that 

v(/i • S - x f 2 ) < |/i| - n and x f 2 = [fi • S]. (5) 

Let {q^/q[ l ' > : i > 0} be the partial quotients of x~ x S_. I 11 [2] (with S = Si, S 2 , ... and 

S_ = J2i>i SiX~ l e x^ 1 F[[x^ 1 ]]) Theorem 1, loc. cit. shows that 

(i) if + {q^l < n < + |gi ,+1 ^| then {qf 1 ,q^) solves (J5]); 

(ii) is an ’n th minimal polynomial’ of S, 0 p. 39]; 

(iii) all ?z th minimal polynomials of S can be expressed in terms of q[^ and q[‘ 11 . 

2.2 Linear Recurring Sequences 

For an infinite sequence S over F, we easily have 5 = 6 F(i) for some if e F[x] 

with 1-01 < \ip\ = d if and only if (<p • S_)i = 0 for i < 0 and xif = [(/?•£]. Now 

(if ■ S)i = ifdSi-d + (pd-iSi-d +1 -f-h (foSi so that S 0 ,..., S^ d , ip and the equation 

Si- d = -( ip d -iSi- d+1 H-h ipoSi)/(p d for i < 0 (6) 

uniquely determine all subsequent terms of S] ip is called a characteristic polynomial of 
the linear recurring sequence S. It is well-known that these polynomials form a (principal) 
ideal Ids of F[x], generated by a minimal polynomial of S. 

The situation is similar for n > 1 and a finite sequence s = so,..., s\- n over F i.e. 
Si G F. If ip G F[x], 1 < d — \ip\ < n and 

Si- d = -(ipd-iSi-d+i H-b ipoSi)/ip d for d + 1 - n < i < 0 

then ip and So,..., si- d uniquely determine S - d ,..., Si_ n and ip is often called a ’character¬ 
istic polynomial’ of s. However, these do not form an ideal of F[x]; the reader may easily 
find examples for which a sum of characteristic polynomials of s is not a characteristic 
polynomial. 


2.3 Annihilating Polynomials and Solutions 

Let n > 1 and s = s 0 ,..., Si_ n be a finite sequence over D i.e. s* G D; s is trivial if 
s = 0,..., 0 = 0 n . The generating function of s is s = Y^i=i-n SiX% e D[x -1 ]. We put 
v = v(s) if s is understood, so that if s is non-trivial then 1 — n < v < 0 . 

Definition 2.1 ( ITU Definition 2.7]) Let p G R, d = \p\ and p = ^2i = oPiX l . If n > 1 
and s = So,..., Si_ n then p is an annihilating polynomial of s, written p G Ann(s), if 

(p ■ s)i = pdSi-d + Pd-iSi-d +1 H-b p 0 Si = 0 for d + 1 - n < i < 0. (7) 

We note that (0) is vacuously satisfied if d = \p\ > n, so that the previous definition is 
equivalent to HB Definition 2.7] and 

{p G R x : | p | > n} C Ann(s) x = {p G R x : (p ■ s)i = 0 for \p\ + 1 — n < i < 0}. 

We prefer ’annihilating polynomial’ to ’characteristic polynomial’ as we do not insist 
that lc(y>) be a unit of D. Further, we may be unable to express si_ n as a linear combi¬ 
nation of So,..., s 2 - n \ e.g. if s = CP -1 ,1 then Si_ n is not a linear combination of 0 ’s; if 
D = Z and s = 2,1 we cannot express 1 as a multiple of 2 in D. 

As in ([5]) we now have: 

Proposition 2.2 For n > 1 and a sequence s = s o, • • •, si_ n over D 

h e Ann(s) x if and only if /i G R x , v(A ■ s - x f 2 ) < |/i| - n and x f 2 = [fi ■ s]. 

We say that f = (/i, f 2 ) G R x x R (or x f 2 /fi if ffi 1 exists) is a solution for s. 

If (p,if) G R 2 is a solution for s and d = \p\ > 0 then |^| = v + d — 1. Also p 
and s 0 ,..., Si_d determine if, since for 1 < i < v + d we have i = (x'fi)i = (p ■ s]j = 
o Fj s i—j where 1 — d < i—j < 0. We include a proof of the following for completeness. 


Proposition 2.3 Let ( p , if) be a solution for s = s 0 ,, s i_ n . If 1 < d = \p\ < n and 
lc(<^) is a unit of D then (i) p~ l G Djjrc -1 ]]; (it) p, s 0 ,..., Si_rf determine s_d, ■ ■ ■, s\- n . 

PROOF. If a — 1 — px~ d /p d G D[x _1 ] then p = p d x d {l — a), a 7 ^ 1 since p 0 and 

\ = ppsp —) = F x ~ i(1 + ° + ° 2 + ■■■)£ Dir 1 ]]. 

Thus v(s — x'fi/p) = v(p ■ s — xip) + v(l /p) < \p\ — n — v{p) = —n, which implies that 
(xif/p)i — Si for n — 1 < i < 0. We know that p and s 0 , ..., Si_d determine if. Hence p 
and s 0 ,..., Si_rf determine s_rf,..., Si_„. □ 

Thus if / is a solution, |/i| < n and lc(/i) is a unit of D then so, ■ ■ ■, Si_ n and fi define a 
linear recurring sequence Sf with S_f — x f 2 /f\ and f\ G Id^j.. 
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2.4 First Examples and the Key Lemma 12.6 

Our first examples include the inductive bases for various results below. 

Example 2.4 (i) The (finite) geometric sequence of length n with common multiple m G 
D x is given by Sq = 1 and Si = m~ l for 1— n < i < — 1. It will be convenient to allow 
m = 0 as well. Then v = 0 and s = m n ~ 1 x 1 ~ n + • • • + 1. If \x— m\ + l— n = 2— n < i < 0 
then ((x — m) ■ s)i = ( x ■ s)i — (m ■ s)i = s;_i — msi = 0, so that x — m G Ann(s) x . As 
x — m is invertible, the corresponding solution is xj[x — m). 

(ii) Let n > 2, a G D x and s = 0 n_ 1 ,a. We have s = ax 1 ~ n , x n = x 1-v G Ann(s) x 
and ax/x n is a solution for s, but cannot express s i_ n as a linear combination of zeroes. 

(in) Let D = Z and s = 2,1- We have the solution (2x — 1, 2) but cannot express s_i 
in terms of sq. 

For s = so, • • •, si-n and arbitrary a G D, t = s, a is the sequence so,..., si_„, a. Given 
a solution g for s we want to construct a solution h for t. We begin with first components. 
It is clear that Ann(f) x C Ann(s) x . Suppose that g\ G Ann(s) x . Then g\ G Anri(f) if 
and only if ( g x ■ t)i = 0 for \gi\ — n < i < 0. If \gf\ + 1 — n < i < 0 and 0 < j < \g\\ then 
1 — n < i — j <0 and so t l _ ] = Sj_j. Hence (g± ■ t)* = (gi ■ s)i — 0 for |(?i | + 1 — n < i < 0 
and gi G Ann(t) if and only if (gi ■ t)\ gi \- n = 0. 

Definition 2.5 ffTIi Definition 2.10], cf. ]?)). Let n > 1, s — s 0 ,..., si_ n and gi G R x . 
For arbitrary a G D and t — s,a the discrepancy of gi and t is A(gp t ) = (gi ■ t)\ gi \- n . 

In general, if g\ G Ann(s) then gi - t = G + A(pi; t ) x^~ n + [gi ■ t] where v(G) < |( 7 i| — n. 
We recall the proof of the following lemma from na as it shows the usefulness of second 
components. Also, the polynomial g 2 hi — g\ h 2 of the proof will reappear later. 

Lemma 2.6 ( JTR Lemma 5.2], cf. Theorem 1]) Let n > 1, s — sq ..., s i_ n and g\ G 
Ann(s) x . Ift = s, a and gi fL Ann(t) then for any hi G Ann(t) x we have \h x \ > n+1 — \gi\. 

Proof. Let A = A (gi;t) 0. We have gi-t = G + Ax^~ n + xg 2 where v(G) < \g± \ — n 
and x g 2 = [gi ■ t\. Also hi ■ t — H + xh 2 where v(Lf) < | h\ \ — n and xh 2 = [h 1 • t\. Put 
(p = g 2 hi — g± h 2 G R. Then 

xip = (xg 2 ) hi - gi (xh 2 ) 

= (gi-t-G-A x^~ n )hi - gi(hi ■ t - H) = -G hi + g x H - A x l5l| ~ n hi 

where v(—Ghi + giH) < |(yfi| — n + |hi| and x<p 7 ^ 0 since D has no zero divisors. Hence 
\gi\ + \hi\ — n — \x(p\ > 1. □ 

For Example I2.4f i). if h is a solution for s then | h 1 1 > 1 since s is non-trivial. For Ex¬ 
ample I2.4f iil with n > 2, (1, 0) is a solution for 0,..., 0 but not for s since a 0, so if 
h is a solution for s then \h\ \ > —v + 1 = n. This can also be proved directly, see [ITl 
Proposition 3.5(c)]. 
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2.5 Linear Complexity and Minimal Solutions 

Next we discuss minimality. Firstly, Arm(s) x 7 ^ 0 since any polynomial of degree n 
annihilates s and the following definition makes sense. 

Definition 2.7 ( ITR Definition 3.1]) Let n > 1 and s = sq. .. ,si_ n . The linear com¬ 
plexity s is 

L n = L(s) = mind/!| : h <G Ann(s) x }. 

We say that f\ is a minimal polynomial (MP) of s if f\ G Ann(s) x and |/i| = L(s). We 
write MP(s) for the set of minimal polynomials of s and say that f 6 R 2 is a minimal 
solution for s if it is a solution for s and fi G MP(s). 

It is important to note that linear complexity and minimality are defined independently 
of how solutions are obtained. Of course, 0/1 is a minimal solution for any sequence of 
zeroes, and L(s) = 0 if and only if s is trivial. For Example I2.4f ii). x n G Ann(s), so 
L„ = n by Lemma 12.61 The function L is a non-decreasing function of n and L n < n. 
If h[ G Ann(s, a) and \hi\ = L(s) then h\ G MP(s, a) since L(s, a) > L(s) = |/q| > L(s, a). 

We repeat the proof of the next result from Em for the convenience of the reader. 

Proposition 2.8 m Corollary 3.24]) Let n > 1, s — s 0 ,..., s i_ n be a sequence over 
F . f be a solution for s and d = gcd(/i, f 2 ). Then (i) f /d is a solution for s; (ii) if f is 
a minimal solution for s then \f±\ < n and d = 1. 

Proof. Let fi ■ s — F + xf 2 where v(F) < \fi\ — n and g — f/d. Then F/d G F[[x -1 ]] 
and gi ■ s — F/d + xg 2 where v(F/d ) = v(F) — \d\ < \fi\ — n — \d\ = |pi| — n. Hence 
gi G Ann(s) and \g± ■ s] = xg 2 i.e. g is a solution for s. (ii) If / is a minimal solution then 
l/il = L n < n and |/i| = L n < l^l = |/ x | — \d\ and hence \d\ = 0. □ 

So if s is a sequence over F, / is any minimal solution for s and S_j = xf 2 / j \, Id^ = fi F[x] 
i.e. /1 is a minimal polynomial for Sf. This justifies our use of the term ’minimal 
polynomial’ of s. The converse of Proposition I2.8f ii) fails: let n > 2 and ,s t = m~ l for 
1 — n < i < 0 where m G F x ; x — m G MP(s) and L n = 1. Now fi = x n G Ann(s), 

f 2 = +-f m n_1 and gcd(/i, f 2 ) = 1 since m n_1 7 ^ 0, but j\ <4 MP(s). 

We will need to single out two kinds of non-trivial sequences: 

Definition 2.9 Let n > 1 and s = so,..., si_ n be a sequence over D. We call s pseudo¬ 
geometric if L n = • • • = Li = 1, and essential if n > 2 and L n > Li > 0. 

Any geometric sequence s is pseudo-geometric since x — S-i/sq G MP(s). A non-trivial 
sequence s = a, 0 n_1 is pseudo-geometric. In general, s is pseudo-geometric if and only if 
s 0 x — s_i G MP(s). Conversely, if n > 2, s 0 is a unit of D and s is pseudo-geometric then 
s is a geometric sequence with common ratio s_i/so- 

Essential sequences were motivated by the need for the integer 

n' = max {/ : L,■ < L n } 

l<j<n 
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to be well-defined; now s' = s o, • ■ ■, si-n' is a well-defined, proper subsequence of s as 
1 < n'< n. 

N.B. For n> 2, the sequence s = 0 n l ,a of Example \2.f\ (ii) is essential since L n = 
n > Li = 0, 1 < n' = n — 1 < n (and moreover s' = CP^ 1 has minimal solution 0/1). On 
the other hand, if s is pseudo-geometric then n' is undefined. 

We can now formally state our subdivision of sequences. 

Proposition 2.10 A sequence over D is either trivial, pseudo-geometric or essential. 

PROOF. We have Lx = • • • = L„ = 0 if and only if s is trivial. Hence if s non-trivial, Lj > 0 
for some i, 1 < i < n. If Lx ^ 0 then Lx = 1 since Lj < i for all i. If 1 = L 2 = • • • = L n 
then s is pseudo-geometric; otherwise if i is the first integer with Lj > 1 then L n > Lj > 1 
i.e. s is essential. Finally, if Lx = 0 then Lj ^ 0 for some i, 2 < i < n as s is non-trivial, 
so L n > Lj > Lx = 0 and s is essential. □ 

Proposition 2.11 If n > 2, s = s 0 , ...,si_ n is pseudo-geometric and // 1 e MP(s) 
satisfies A(/xx;s, s_ n ) ^ 0 then s,s- n is essential. 

PROOF. By Lemma [2.61 L n+ x > n + 1 — L„ = n > 2 > 1 = Lx. □ 

The first state diagram of the Introduction illustrates this transition on adding a term. 


3 Minimal Solutions via Partial Quotients 

Here we revisit the first part of [9j Theorem 1], Let S be an infinite sequence over F, S_ 
its generating function and n > 1. Our goal is to show that a certain partial quotient of 
S_ (depending on n) is a minimal solution for S\n — S 0 ,... , Si_ n . In particular, we relate 
L(5'|n) to the degrees of the partial quotients of S_. 

We recall the construction of the partial quotients of S_, their basic identities and 
properties. We work through m Example 1], Then we discuss geometric sequences and 
0 Tl ^ 1 ,5'i_ n , where Sx-n £ P x and n > 2. These form our inductive basis for the main 
Theorem 13.71 When S\n is essential, we determine maxx<j <n {Lj < L n } and prove an 
identity for any / 6 R 2 . 

3.1 Continued Fractions 

We use the formulation of continued fractions in F[[x -1 ]] from [9]; in particular, we also 
use Pol(L) = Y7i=o L i xl G R for L G F((x 1 )). It is well known that S_ has the unique 
continued fraction expansion 0 + x/(a\ + l/(a 2 + • • •)) = [0, a±, a 2 ,...] where, if aj G R 
exists, then |aj| > 1. The aj are obtained using division in the field F((a: _1 )) as follows: 
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a 0 <- 0; b 0 4 - x 1 S ; q { 1] «- (0,1); g (0) «- (1, 0); 


i i — 0; 

while [hi 7^ 0) 


I" a i+ 1 Pol(6j x ); q (i+1) a i+ iq w + g ( * x) ; 

6j + i fej 1 — a i+1 ; i •<— i + 1; J 

The partial quotients of S_ are { q^ 6 R 2 : i > 0}. Put | q[ ^| = 0 and if bi = 0 (i.e. 
a i+ i and f/* + P do not exist) put |a i+ i| = |qj* +1 ' l | = oo. In the following well-known result, 
Part (iv) on numerators is probably well-known, but does not appear in [9|. 

Theorem 3.1 If i > 1, q — q^> exists and q' = gh-i) then \qi\ = )T)' =1 \aj\ > 1. In 

particular, if i> 0 then 1 — |gf +1 | < 0 . 

If i > 0 and q exists then 

(i) Q 2 d'l ~ qi 92 = (- 1 )*' 1 and gcd(gi, q[) = gcd (q 2 , q' 2 ) = 1; 

(a) 

_ x (q 2 + bi q' 2 ) 

~ qi + biq[ 

(Hi) v(q\ ■ S_ — xq 2 ) = 1 — |gf +1 ^| so that Pol(q! ■ S_ — xq 2 ) = 0; 

(iv) |Oi| = 1 - V, |< 72 1 = Y!j= 2 l a il = kil - l a i| ifi > 1 and Pol(gi • S) = [g x ■ S] = xq 2 . 

If i > 1 is the first index for which bi = 0 then 

(v) S_ = x q 2 /qi and v(q 2 ■ S_ — xq 2 ) = —oo = 1 — |gf +1 |; 

(vi) Pol(gi • S_) — xq 2 and \q 2 \ = v + |gi| - 1 < |gi|. 

PROOF, (i)-(iii) Use induction and properties of the exponential valuation as in [9] . (iv) 
For i > 1, |g 2 | = Y2]= 2 \ a j\ * s an eas Y induction. We have |ai| = |Pol(x/iS)| = 1 — v so 

that || = |<?i| - |Oi| = |?i| + v - 1 and 


v+ki| I92I 

Pol(q , i ■ S) = Pol(qi ■ S_~xq 2 ) + ^ (x q 2 )j x 3 = x ^ ( q 2 )kX l = xq 2 

j =0 k =—1 

by (iii). (v) If i > 1, 6 ,;_i 7 ^ 0 and bi = 0 then 

Q = Q + (Kl 1 - ai)q =q- a,4 + b~\q = g ( ^ 2) + bf\q 


and rearranging gives 

xq 2 = x(q! 2 + b i - 1 q$~ 2) ) _ g 
Qi q[ + b i _ 1 q{ 2 ~ 2) ~ 

by (ii). (vi) Immediate. 


□ 
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Next we define a partition of N x . Let no = 1 and for i > 1, define r rii = nfiS ) by 

ni = k!* _1) | + 

From Theorem 13.11 we have 1 = no < n\ < U 2 < ■ ■ ■ and {[ni, nj+i) : i > 0} is a 
partition of N x (except that [no, ni) = 0 if ni = 1) for if i G N is largest such that n* < n 
then n G [rij, nj+i) and i is clearly unique. Thus if S_ = 0, ri\ = oo and for all n > 1, 0/1 
is a minimal solution for S\n. 

The next example is m Example, p. 21]. Here F = F 2 [a;]/(a; 4 + a: + l) and a generates 
F x . For the table of F x as polynomials in a , see HU p. 85], 

Example 3.2 Let S_ = x(a 5 x 2 + a 2 x + a 10 )/(x 3 + a 6 x 2 + a 3 x + a 13 ) G ¥(x). Direct cal¬ 
culation gives 


i 

k 

^Z+l 

^Z+l ^ ^2+1 

0 

x~ 3 S 

a 10 x + a 14 

{a^x + a 10 ) / (a 5 x 2 + a 2 x + a 10 ) 

1 

(cfix + a 10 ) / (ofix 2 + a 2 x + a 10 ) 

a u x + a 5 

a 5 /{a G x + a 10 ) 

2 

a 5 / ( a 6 x + a 10 ) 

ax + a 5 

0 


i 

gb“b 

gW 

gb+i) — ai+1 gb) -(- gb b 

0 

(0,1) 

(1,0) 

(< a 10 x + a 14 ,1) 

1 

(1,0) 

(a 10 a; + a 14 ,1) 

( a 9 x 2 + a 6 x + a, a 14 a; + a 5 ) 


and q ® is 

(a 10 a: 3 + a 13 x 2 + a 3 x + a 8 , x 2 + a 12 x + a 5 ) = a 10 (x 3 + a 6 x 2 + a 3 x + a 13 , ot'x 2 + ot 2 x + a 10 ). 

Clearing denominators, this is the extended Euclidean algorithm: a* + i is the quotient and 
b i+ i is the remainder. We have = 0, |g^| = 1, |g[ 2 ^| = 2 and |g[ 3) | = 3 and 
|c/i 4) | = oo, so that the partition o/N x defined by S is [1,3), [3,5), [5, oo). 

Our inductive proof of the first part of 0 Theorem 1] depends on characterising 
solutions for S\n in terms of solutions for S, and is proved using Proposition 12.21 

Lemma 3.3 Let D = F be a field. If f\ G R x , |/i| < n and s = S\n then [f\ ■ S] = 
Ifi ■ s] = x f 2 say. Further f is a solution for s if and only if v(/i • S_ — x / 2 ) < |/i | — n. 

Proof. We have v(5 — s) < —n and [fi ■ (S_ — s)] = 0, so [fi ■ S] = [fi ■ (s + (S_ — s))] = 
[fi ■ §\ + [fi ■ (S - s)] = [fi ■ s] and 

fi-S-x\fi-S] = fi-(s + (S-s))-x{f 1 -s} = f!-s-x[fi-s} +fifiS-s). 

Hence v(/i ■ S_ — x[f\ ■ S]) < max{v(/i • s — x[fi ■ s]), |/i| — n} and if fi G Ann(s) then 
v (/i • S_ — [f\ ■ A]) < |/] | — n by Proposition 12.21 The converse is proved similarly, for 
fi-s~x[f 1 .s\ = f 1 .S~x[f 1 .S]-f 1 .{S-s). □ 
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3.2 Geometric Sequences 

If So 7 ^ 0, Theorem 13.11 implies that n\ = |ai| = 1 — v = 1. 

Proposition 3.4 Let S be an infinite sequence over a field F such that So 0, q — q 1 ' 1 ' 1 
and q' = (1, 0). The following are equivalent 

(i) ne [ni, n 2 ); 

(ii) q' is not a solution for s = S\n and q is a minimal solution for s. 

PROOF, (i) =>■ (ii). Firstly q' is not a solution for s since S 0 0. Secondly, 1^1 < n and 

v (*?i • S_ - xq 2 ) — 1 - \q^\ — 1 - (n 2 - 1 ) = 2 - n 2 < 1 - n. 

By Lemma [3.31 q is a solution for s. As |gi| = 1, it is a minimal solution. 

(ii) =>■ (i). If n ^ [ 711 , 712 ) then n 6 [ 712 , 713 ) so v(qi-S_—xq 2 ) = 1 —1(^| = 2 — 7i 2 > 1—77 
and q\ ^ Ann(s). In particular, q± ^ MP(s). □ 

We can say more. 

Proposition 3.5 Let S 0 7 ^ 0 and r = S'_ 1 /S' 0 . Then ai — (x — r)/S 0 and n 2 > 3. 

Proof. Write T_ = S_/S 0 = 1 +rx~ 1 + - ■ ■ ; if this is a geometric series then T_ = x/(x — r), 
bo = x~ 2 S_ — So/(x — r), ai — PoI^q 1 ) = (x — r)/S 0 and n 2 = 00 . 

Otherwise we have T — 1 + rx~ l + r 2 x~ 2 + • • • + r n2 ~ 2 x 2 ~ n2 + T x _ n2 + • • • where Ti_ nz 7 ^ 
rT 2 _ n2 and 2 — n 2 < —1. Now T = x/(x — r) + U_ for some LJ_ with v(t/) < 1 — n 2 < — 2 
and 

t , 0 -1 p/p - 1 ^ t i-i , -i rr 1 + (x - r)^- 1 !/ 1-y 

bo So = x SI Sq — x T_={x — r) +x U_ — -=- 

a: — r x — r 

say, where v(V) < -2. Thus for all i > 1 we have v(W) < -2 i < -2, 

60 1 = ^ 0 _ 1 (^-^)/(l-K) = So\x-r)(l + V + V 2 + ■■■) 
and a\ = Po1(5q *) = (x — r)/So . □ 


3.3 Essential Sequences and the General Case 

If So = 0, Theorem 13.11 implies that 77i = |ai| = 1 — v > 1. 

Proposition 3.6 Let S be an infinite sequence over a field F such that S_fi0, So = 0, 
q = q^ and q' = (1, 0). The following are equivalent 

(i) 7i G [ 711 , 712 ); 

(ii) q' is not a solution for s = S\n and q is a minimal solution for s. 
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Proof, (i) => (ii). Let n G [ 711 , 712 ). We have |gi| = 1 — v = ni<n and 
y(qi ■ S-xq 2 ) = 1 - |g| 2) | = 1 + |gi| -n 2 < \qi\ -n 

from Theorem 18.11 so that q is a solution for s by Lemma 18.81 As nj - 1 > 1 and 

1 G Ann(5'|ni — 1) \ Ann(5'|ni), Lemma 12.61 implies that 1 — v = \qi\ > L n > L ni > 

711 — L ni _i = Tii = 1 — v so L„ = 1 - v and q is a minimal solution for s. 

(ii) =>• (i). If n 0 [ 771 , 712 ) then either (a) n G [l,77i), v(q[ • & — xq 2 ) = 1 — |gi| = 
1 — Til < —n and q[ = 1 G Ann(s) or (b) n G [ 712 , 77 . 3 ) and as in the proof of Proposition 
13.41 v(gi • S_ — xq 2 ) = 1 — |g[ 2) | = 2 — 77 2 > 1 — n, so q\ ^ Arm(s) and q\ fL MP(s). □ 

We have now treated the case n G [ 771 , n 2 ). Now for the general case. 

Theorem 3.7 (Cf. f9, Theorem 1], [El Theorem 4]) Let S ^ 0, j > 1 and q = q^>, 

q' = gh-i). The following are equivalent: 

(i) 7i G [ni,n i+ i); 

(ii) q' is not a solution for s = S\n and q is a minimal solution for s. 

Proof. For 7 = 1 the result follows from Propositions 13.41 and 13.61 Suppose inductively 
that 7 > 2 and that the result is true for 7 — 1. 

(i) (ii). Let n G [ 77 , 4 , 77 * 44 ). Then 0 < |g[| = 77 * — |gi| < n — \qi\ i.e. |gi| < n and 
Lemma [3.31 applies. If n < 77* +1 then by Theorem 13.11 

v(qi ■ S_-xq 2 ) = 1 - |gf +1) | = 1 -77*41 + |<?l| < ki| -77 

so gi G Ann(s). Likewise \q[\ < |gi| < n and Lemma [3731 applies to q[: if 77 * < n then 

v(g[ • S-xq 2 ) = 1 - \qi\ = 1 - 77 * + \q[\ > 1 + \q[\ - n 

i.e. q[ ^ Ann(s). We next show that gi G MP(S'|t 7 *). Since gi G Ann(s) we have 
gi G Ann(S'|n*), so |gi| > L ni . We have 77* > 77 2 > 77 i > 1 i.e. 77* — 1 > 1. The inductive 
hypothesis and Lemma [3.31 imply that q[ G Ann (iSI 77 * — 1) \ Ann(S'| 77 *). So Lemma [2.61 
implies that |gi| > L ni > 77* — \q[\ = |gi| and hence |gi| = L ni . Now let 77* + 1 < 77 < 77*41. 
We know that q± G Ann(s), so L n < |gi| = L ni . But n > 77* implies that L n > L***, so 
L** = L n . and gi G MP(s). 

(ii) =^- (i). If 77 ^ [ 77 *, 77 * 41 ) either (a) n G [77*_i,77*) and q[ G Ann(s) or (b) n G 

[ 77 * 41 , 77 * 42 ) and so gi 4 L Ann(s) by the first part, and in particular q x 4 L MP(s). □ 

It follows that for Example 13.21 we have Li = L 2 = 1, L 3 = L 4 = 2 and L 5 = Lq = 3. 
We list some simple consequences of Theorem 13.71 

Proposition 3.8 Let n > 1, s — S\n be non-trivial. Then 
(i) either s is geometric or essential; 

(n) L n . = |afc| and we can obtain L n . without computing partial quotients; 

(Hi) 7 / 77*44 < 00 then on the interval [ 774 , 77 * 41 ), L n . = |gj^| appears n i+ \ — ? 7 * times. 
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Proof. We prove (i) only. As s is non-trivial, n [1, ri\) i.e. n £ [rii, rii+ 1 ) for some i > 1 
and ri\ = |g x | = |a x | = 1 — v. If i > 2 then L n = L n . > L n2 > L ni > L x i.e. s is essential. 
Now suppose that i — 1 i.e. n £ [n x , n 2 ). If v = 0 then |g x | = 1, L n = L ni = L x = 1 and s 
is geometric. If v < 0, ri\ = |g x | > 1 and L n = L n . > L ni = 1 — v > 2 and L ni > L x = 0 
i.e. s is essential. □ 

We note that in the previous proposition, if n £ [n i, n 2 ) and fe x = 0 (i.e. n 2 = oo) then 
S'ln is geometric; otherwise for k £ [n.j,nj +x ) and i > 2, S'l/c will be essential. The second 
state diagram of the Introduction illustrates this behaviour. 

Corollary 3.9 (Cf. J7| Theorem 1]) Let n > 1 and s = S'|n be non-trivial. Then 

(i) either L n+X = L n or L n+ i = n + 1 - L n > L„; 

(ii) if n £ [rii,ni- |_i) and gf ^ Ann(5'|n + 1) then L n+ i = max{L n ,n + 1 — L n }. 

PROOF. We know that n £ [ni,Ui +x ) for some i > 1. (i) From Theorem 13.71 if n + 1 G 
[n;,n i+ i) then L n+ i = L Ui = L n . Otherwise n + 1 = n , ;+1 = L n . + L n . +1 = L n + L n+ i and 
L n+ i = L n . +1 > L„. = L n . (ii) If L n > n + 1 - L„ then L n+ i = L n by Part (i). Suppose 
now that L n < n + 1 — L n . Since n + 1 — L n < L n+ i by Lemma 12.61 L n < L n+ i and by 
Part (i) we have L n+ i = n + 1 — L n . We conclude that L n+ i = max{L n , n + 1 — L n }. □ 

If S\n is essential, the integer ma ~X-i<j <n {j : Lj < L n } equals nj — 1: 

Corollary 3.10 If L n > L x > 0, n £ [rii,ni +x ) ; s = 5|n, g = gW } cf = gh-b and 
n' — n t — 1 then 

(i) cf is a minimal solution for s' = S\n'; 

(ii) n' = ma xi<j <n {j ■ L j < L n } and L n + L n / = n' + 1; 

(in) [g x • s] = xg 2 and [q[ ■ s] = xq' 2 . 

Proof. We know that n £ [rq,nj + i) for some i > 1 and L n . = L n > L x > 0. Hence 
rf > 1 and s' is well-dehned. As n' £ g( G MP(s') by Theorem 13.71 (ii) We have 

L n ' = L ni _i < L ni = L n and so n' = maxi< i<n {j : Lj < L„}. Also L n + L n / = |g| + |g'| = 
Hi — n' + 1. (iii) We have [g x -s] = [g x • *S] = xg 2 by Lemma [3731 since |g x | = L ni <ni<n. 
Likewise |g(| < |g x | < n and [q[ ■ s] = [q[ ■ S\ = xq' 2 . □ 

We conclude with a consequence of Theorem 13.11 and Corollary 13.101 

Proposition 3.11 Suppose that L n > L x > 0, n G [iii,n i+x ), g = gW and cf = q^ l ~ l \ If 
f £ R 2 , m = f 2 g x - fi g 2 and ml = f 2 q[- /, g^ then 

( 1)* —1 f\ — m! qi — m q[. 

PROOF. As cf is well-dehned from Corollarv l3.101 we have ml g x — m q[ = (/ 2 g) — f\ q 2 ) g x — 
(/2 gi - fi qf) q'l = /i(g 2 g( - g 2 'gi), which is (-1 )* -1 /i by Theorem O □ 
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4 An Inductive Construction of Minimal Solutions 


In this section, we work with arbitrary finite sequences over D. Given n> 1, we construct 
a minimal solution // for any s — s 0 , ■ ■ ■, Si_ n i.e. H\ G Ann(s) and |/xi| = L n = L say. 
If L < n and lc(//i) is a unit of D, we can ’generate’ s_l, ..., si_ n using s 0 ,..., s x _l and 
/i]. For Examples I2.4f ii). (iii) we will see that the construction returns Hi = x n and 
Hi — 2x — 1 respectively; in each case, we cannot generate Si_ n using s 0 ,..., s 2 _ n and Hi- 
We simplify [TT] by appealing to Proposition 12.101 and recalling two lemmas from [ITj. 
In this way we construct a new minimal solution when the current one fails. The proof 
of each lemma consists of (i) verifying that we have a new solution and (ii) applying 
Lemma 12.61 to deduce minimality. For a pseudo-geometric sequence, it suffices to consider 
n + 1 — 2L = 77 — 1 >0 only and the proof is elementary. However, for an essential 
sequence, we require both a current minimal solution and a solution for s 0 ,..., S\_ n ' 
where 1 < n' < n. We encode each of these solutions as a ’triple’. 

The resulting Algorithm 14.121 is identical to [IT, Algorithm 4.6], except that we com¬ 
pute /i -2 rather that xfii- We can also suppress second components and in this way 
compute minimal polynomials only, cf. [7]. We include a normalised version to compute 
a monic Hi when D is a field. We also give a new analysis of Algorithm 14.121 

Section [25 defines the scalar V s G D x and proves Identities (j2J), © of the Introduc¬ 
tion, see Propositions 14.20114.221 These identities are integral to the rest of the paper. 

4.1 Pseudo-Geometric Sequences 

The following integer will play an important role for all finite sequences. 

Definition 4.1 For n > 1 and s = sq, ..., s x _ n we put e s = 77 + 1 — 2L n e Z. 

Lemma 4.2 m Theorems 3.8, 4-5]) Let k > 1, r = s 0 ,...,si_fc, so 7 ^ 0 and h be a 
minimal solution for r with L^ = |/ii| = 1 . If s = r,S-k and A x = A(hi',s) 7 ^ 0 then 
v = sox er H ~ Ai(l, 0) is a minimal solution for s; in fact |zq| = max{l, 1 + e r }. 

We apply this as follows. If 77 = 1 and so 7 ^ 0 then /j = (x, sq), a pseudo-geometric 
sequence and e so = 0. For n = 2, r = s 0 and s = r, s_ 1; if A x = s_ x 7 ^ 0 then v = 
s 0 x°(x, s 0 ) — s_i(l, 0) = (s 0 a: — s_i, Sq) since e r = 0; s is also pseudo-geometric. Let 77 > 3 
and s = r, s x _ n . If L n _ x = • • • = L x , ^ is a minimal solution for r and A x = A (ht, s ) 7^ 0, 
then v = soX n ~ 2 H — A x (l, 0) since e r — n — 2. Now L n — \v x \ — n — 1 > L x = 1; the new 
sequence s is essential. This is an explicit version of Proposition 12.111 

Examples 4.3 Let a,b E D x . (i) Put r = a and s — r,b. We begin with /j = (x, a ) and 
A x = A (h; s) = b 7 ^ 0. Lemma shows that a(x, a) — 6(1, 0) = (ax — 6 , a 2 ) is a minimal 
solution for s. If a is a unit of D and m = 6/a, we have s = a,ma i.e. Example \2.^( i) 
and since ax — 6 G MP(s), we have x — m G MP(s). 

(ii) Now let k = 2, r = a, 6 and s — r,c. We know that (ax — b,a 2 ) is a minimal 
solution for r and 

A x = A(ax — 6; s) = ((ax — b) ■ (cx~ 2 + bx + a) ) x _ 2 = ac — 6 2 . 
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Lemma \J7^ implies that if Ai 7 ^ 0, ax(ax — 6 , a 2 ) — Ai(l,0) = ( a 2 x 2 — abx — Ai,a 3 x) zs 
a minimal solution for s. Taking D = Z, a = b = 1 and c = 2, we see that x 2 — x — 1 G 
MP(1,1, 2), as expected. Further, if f\ = x 3 f/zen f 2 = x 2 + x + 2 and gcd(/i, / 2 ) = 1. 
T/zzzs t/ze converse of Proposition HOI fails for essential sequences too. 

(Hi) Let the common multiple m in Example \2.f\ (i) be zero, so that r = a, where 
k >2 and // = (x,a) i.e. e r — k — 1. If s = r,b where b G D x t/zen Ai = A (a;; s) = 6 7 ^ 0, 
/zence ax k ^{x, a) — 6(1, 0) = (aa^ — 6 , a 2 x k ~ 1 ) is a minimal solution for s and L fc+1 = k > 
1 = L*.. 


4.2 Essential Sequences or, a Tale of Two Triples 

Next we recall a lemma which constructs a minimal solution for an essential sequence 
when the current one fails. As this is more involved, we encode the data as a ’triple’ 
consisting of a strictly positive integer, a minimal polynomial and an element of D x . We 
also require that our two triples are linked by linear complexity. Thus given a pair of 
linked triples for r and a G D, we construct a pair of linked triples for s = r, a. 

Definition 4.4 Let k > 1, r — s 0) ■ ■ ■, si -k be essential and s_ k G D. A linked pair of 
triples (for r) consists ofT r = (k, jx 1 , Ai), T' r = (k', , A() e N x x R x D x such that 

(i) s — r, s- k , Ai = A(/ii; s), gi 6 MP(r) \ Ann(s); 

(ii) k' = maxi<j<fc{j : Lj < Lfc} < k so that r' — sq, , s 1 -k 1 is a well-defined, proper 

subsequence of r; 

(Hi) s' = r', S-k', A( = A(/z' 1 ; s'), p! x G MP(r') \ Ann(s') and + L & = k' + 1. 

Example 4.5 Leta,beD x . 

(i) Let r = a, 0 fc_1 with k >2 (a pseudo-geometric sequence with common multiple 0 ), 
s = r,b and v = ( ax k — 6, a 2 x k ~ 1 ) be a minimal solution for s as in Example \4 ■ iii). Let 
T’ s = ( k,x,b ) and suppose that A x = A(zq; s, c) 7 ^ 0 for some c G D. Then T s — (k + 
1, z/i, Ai) and T' are linked: for i < k we have L, : = 1, so (fc +1)' = k and Lfc + i + L(*, + i y = 
k + 1 = (A; + 1)' + 1. 

(ii) Let k >2 and r = O^” 1 , a where k >2, so v = 1 — k. We claim that T r = (k, x k , b) 
and T' r — (k — 1,1, a) are linked. We have A(l;0 fc ' _1 ) = a 7 ^ 0, giving the triple T' r = 
(—v, l,a) i.e. kl = —v = k — 1 and L k > = 0. We also know that ( x k ,a ) is a minimal 
solution for r. Let s = r,b where b 7 ^ 0. Then 

A(x k ] s ) = ( x k ■ s) k _ k = s_ k = b 7 ^ 0. 

giving the triple T r = ( k , x k ,b). Furthermore T r , are linked since L k + L*,/ = A: = k' + 1. 

Remark 4.6 In m we used T r = (l,x, Ai), T' r = (0,1, so), linked by L 0 = 0 when 
r = s 0 7 ^ 0. In this paper, we treat pseudo-geometric sequences separately and the proper 
subsequence r' always has length n' > 1 , further simplifying the theory developed in m- 

We now combine several results from HU to construct a linked triple for s = r, s_ k 
from a linked triple for r. 
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Lemma 4.7 m Proposition 3.11, Theorem 3.13, Proposition 4-4]) Let k > 2, r = 
s 0 ,...,si_fc be essential and s = r,s_k- Suppose that T r = (k,pi, Af), T' r = [k!, p[, A' x ) 
are linked triples for r. If 

( A[ p — x~ e p / if e = e r < 0 

”= { 

{ A[ x +e p — Ai pi otherwise 

then \vf\ = maxjLfc, L+ e r } = max{Lfc, k + 1 — L^}, /q G MP(s) and x z / 2 = [v\ ■ t\ i.e. 
v is a minimal solution for s. Further, if 

(a) A = A[yi, s , s_fc_x) 7 ^ 0 and T s — (k + 1, iq, A); 

(b) T' s = T' r if e r < 0 and T' s = [k,p 1 , A x ) if e r > 1; 

then T S ,T' S are linked. Finally e s = 1 + e r if e r < 0, and e s = 1 — e r otherwise. 

PROOF. We show only that T S ,T' S are linked. (The remaining item on updating e s is a 
simple verihcation.) Let = ma '■ L j < L^} and pi G MP(r). If e r < 0 then 
L fc+1 + L( fc+1 )/ = L fc + L fc / = k' + 1 = (fc + l)' + l and (fc + 1)' = k' = = m^+i. Otherwise 

L fc+1 + L( fc+1 y = {k + 1 — L fc ) + L fc = k + 1 = (k + 1)' + 1 and (k + 1)' = k = uik+i since 

Lfc < k + 1 — Lfc = Lfc +1 . □ 


Example 4.8 (%) For Example \4-5\ ii), Theorem \4- 7\ yields = ax k — bx k ~ l from T r = 
( k, x k , b), T' r = {k — 1,1, a) which are linked as we have seen, (ii) In Lemma \4. HI L ]. = 1, 
ui = sox k ~ 1 pi — Ai and L^+i — k. If A = A(up, s) / 0 we have T s — (k + l,vi,A), 
T' s = (k , pi, Ai) and Lfc + i + L(fc + i)/ = k + 1 = {k + 1)' + 1 so that T s , T s / are linked. 


4.3 The Inductive Theorem and the Corresponding Algorithm 

The elementary case n — 1, Example I4.5l ii). Lemmas 14.21 and 14.71 now yield 

Theorem 4.9 m Theorems 3.13, 4-5]) For n > 1 and any sequence s = s 0 ,..., s i_ n 
over D, we can construct a minimal solution p for s. 

PROOF. We induct on n. For n — 1, p — (1,0) is minimal if s is trivial and otherwise 
p = (x, so) is- Let n = 2, r = sq and s = r, s_i. If Ai = A (pp s) = 0 then p is as required. 
Otherwise Ai 7 ^ 0 and s is non-trivial, so s is either pseudo-geometric or essential. In the 
first case, r is also pseudo-geometric and we can apply Lemma 14.21 to r and s_i: we take 
p = (sox — s_i, Sq). For the second case, r is trivial and s_i ^ 0 so we take p = (x 2 , s_i). 
Morover if we put n' = 1 then L ra + L n i = 2 = n' + 1. Hence if s is non-trivial and 
Ax = A(x 2 ; s, s_ 2 ) 7 ^ 0, we have linked triples T s = ( 2 , x 2 , Ax), T' s = (1, 1 , s_x). 

Now let n > 3, s = r, Sx_ n and p be our solution for r with linked triples T r , T' r if both 
r 7 ^ 0 and Ax = A (pp s) 7 ^ 0. Thus s is non-trivial; if s is pseudo-geometric, so is r and 
we apply Lemma [4.21 to r and Sx_ n . Otherwise s is essential. If r — 0 then sx_ n 7 ^ 0. Put 
p = (x n , Sx_ n ). Now n' = n — 1, L n + L n / = n = n' + 1 and T' s = (n', 1, sx_ n ) is a triple. 
For r 7 ^ 0, the inductive hypothesis and Lemma [4771 apply to r, Sx_ n , p and linked triples 
T r , T' r to yield a new p, and a linked T s , T' s if s is non-trivial and A(pp s, s_ n ) 7 ^ 0. □ 
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Remark 4.10 The proof of Theorem \f.E\ for minimal polynomials only does not require 
the fact that xp 2 = \p\ ■ s]. 

Corollary 4.11 (JfQ Theorem 2], cf. Corollary \3.<A) If s = s 0 ,...,s x _ n is non-trivial 
and p\ Ann(s, s_ n ) then L n+ i = max{L n , n + 1 — L n }. 

Proof. As pi Ann(s, s_ n ), Theorem 14.91 implies that iq G MP(s,s_ n ) i.e. L n+1 = \u^\. 
If e s < 0, L n+ i = L n > n + 1 — L n . Otherwise e s > 0 and L n+ i = n + 1 — L n > L n . □ 

Next we derive the algorithm which follows from the constructive proof of Theorem 14.91 
The constructions in Example I4.5l ii) and Lemma 14.21 bear some resemblance to Lemma 
14.71 and pseudo-geometric sequences often become essential. Thus it is reasonable to try 
to fit these two cases into the format of Lemma 14.71 and to iterate. 

First we rewrite Lemma [4.71 algorithmically using p for the current solution, the vari¬ 
able A', updates for e and v for the new solution. Here 2 < k < n — 1 and s — s 0 , ■ ■ ■, s \-k 
is essential: 

Lemma 14.71 (restated) 

Ai A(pi; So,, S-fc); 

if Ai 7 ^ 0 then if e < 0 then \u (A) p — Ai x~ e p'\ e <— 1 + e; j 

else |V ■<— A[ x e p — A x (//I, A() (/x, A x ); e 1 — e; J 

(i) We observe how Lemma 14.71 (restated) reduces when k = 0 if we start with /j = 

(1, 0), e = 1, /i' = (0, —1) and A) = 1 : we have A(/^ x ; so) = so and thus if s 0 ^0 we have 
v = A( x e fi — A x p! = x(l, 0) — so(0, —1) = (x, so) and the case e < 0 does not arise. We 

have the correct result when s 0 0; now e = 0, v' — (1, 0) and A) = s 0 . 

(ii) Now let k — 2 and put p — is, p! — v'. We have A x = A (pp s ) = s_ x . If s_ x ^ 0 
and e = 0 i.e. s 0 ^ 0 we have v = A \p — A x x~ e p r = So(x, s 0 ) — s_ x (l, 0) = ( s 0 x — s_ x , Sq) 
as desired. But if s 0 = 0 then s is essential and p = (1, 0), giving v = A\ x e p — A ip' — 
t 2 (l,0) — s_i(0, — 1) = (t 2 ,s_i), v' = (1,0) and A) = s_ x , provided e = 2. So Lemma 
14.71 (restated) behaves correctly when k — 2, provided e = 2 if so = 0. We conclude that 
when A x ^ 0, Lemma 14.71 (restated) applies if we initialise as in case (i) and e 1 + e 
if s 0 = 0. Moreover s, s _ 2 will be essential, so that Lemma l4~7l (restated! can be reapplied. 

Now replace v by p and z/' by p' throughout. This requires a temporary variable T 
to avoid overwriting p when e > 0. We factor out incrementing e, giving last statement 
e •<— 1 + e. 

If A x = 0 then p\ remains unchanged and e 1 + e since k + 2 — 2 \pi\ = 1 + e. We 
can thus place e <— 1 + e at the bottom of the loop, independently of A x , as in: 

Algorithm 4.12 (IT1, Algorithm f. 6 ], cf. (7) Algorithm 1]) 

Input: The n > 1 values of a sequence s = s 0 ,, s x _ n over D. 

Output: A minimal solution p for s. 
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f/i (1, 0); pi G- (0, —1),' A[ •<— 1; e G- 1; 


for i 4 — 0 to 1 — 77 . do 


|" Ai A (pp, sq, , s^; 

if Ai 7 ^ 0 then if e < 0 then p G- A' x p — Ai x~ e p 

else | T G- p; p G- A' x x e p — A x pi ; 
(/I A^) G- (T, Ai); e 4 - e; J 


e 1 + e; J 

return /i.J 

Note that after so = 0 we have e > 2. We verify the remaining cases: 

(iii) s trivial; Algorithm 14. 121 gives p = (1,0) as it should; 

(iv) s = 0 fc- 1 ,Si_fc where Si_a, G F x ; here p = (1,0), e = [k — 1) + 1 — 2\pf\ = k 
and A(/xp s 0 ,..., si_ fc ) = S]_ fc , Algorithm 14.121 gives /j, = A[x e p — A x /j' = ;r fc (l, 0) — 
Si_fc(0, —1) = (x fc ,Si_fc), p' = (1,0) and A' x = si_ fc . This agrees with Example 14. 51 Ii) and 
moreover Lemma [4.71 can be reapplied. 

(v) k > 2 and so> •• • , si-fc is pseudo-geometric; here p = (sqx — s_i,Sq), p' = (1,0) 
from (ii) above and e = k + 1 — 2\pi\ = k — 1. If Ai = A(pp, so,..., S-k ) ^ 0 then Algo¬ 
rithm 14.121 gives A[ x e p — A\ p' — A[ x k ~ l p — Ai p' and /i = //, A' x = Ai, which agrees 
with Lemma 14.21 and Lemma 14.71 can be reapplied. Finally, if Ai = 0 then p is unchanged. 


We conclude that Algorithm 14.121 computes a minimal solution p for s. Note that 
(i) we may suppress second components and compute /j i only as in [7]; (ii) Algorithm 
14.121 is identical to HD Algorithm 4.6] except that //' G- (0, —a;) has been replaced by 
fi' G- (0, —1), so that Algorithm 4.6, Zoc. cit. computes x/x 2 instead of p 2 . 


Remark 4.13 (Initialisation) In J2 Section 7.3], (a®,a/ 0 )) = (1,1) and (r^, 7 ^) = 
(1, 0). This corresponds to the fact that 1 + Six + • • • + s n a: n G F[x] is used in the key 
equation fl\, Equation 7.302]. Thus if s = 0 n then 1 = 1/1 obtains in flf, whereas 0 = 0/1 
obtains in our approach. 

Our initialisation pi = (0,-1) was chosen to yield the inductive bases of Theorem 
IE In E’ we have the initialisation 'B(D) — 1 which corresponds to p\ = 1. Let 
Algorithm \ f, 12 \ denote Algorithm \4.12\ using the initialisation p[ — 1. The reader may 
easily check that the first iteration of Algorithm \ f.!2( (with A\ = s v 7 ^ 0) produces 
pi = x n — s v G MP(s). As Lemmas \f.2\ \f. 9| apply to any pi G MP(s), Theorem \4-9\ and 
hence Algorithm\4-12\' also produces a minimal polynomial on subsequent iterations. 
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Tabic 1: Algorithm 14.121 for Example I4.3f ii) 


s 

A, 

a; 

e s 

p 

p! 


— 

1 

1 

( 1 . 0 ) 

( 0 .- 1 ) 

a 

a 

1 

1 

(. x , a) 

( 1 . 0 ) 

a, b 

b 

1 

0 

(ax — b , a 2 ) 

( 1 . 0 ) 

a, b, c 

ac — b 2 

a 

1 

ax(ax — b , a 2 ) — Ai(l, 0) 

(ax — b , a 2 ). 


Proposition 4.14 (Normalised Algorithm 14.121) //D = F is a field, p = Ai/A[ 
and p of Algorithm\£T%is updated via 

f p — p x~ e p! if e < 0 
^ 1 x +e p — p p! otherwise 


then Algorithm \f,13\ produces a minimal solution p for s with pi monic. 


PROOF. It suffices to show that the updating is well-defined and p\ is monic. Firstly, 
A[ = 1 initially and A( is either unchanged or replaced by Ai 7 ^ 0. Thus p is well-defined. 
Secondly, p\ is monic for the base cases. Suppose that s — Sq, , s i_ n is essential and 
e < 0. Then lc (pi — p x~ e p'fi) = lc(^i) since —e + L n / = 2L n — n — 1 + L„/ = L n + n' — n < L n 
as L n + L n , = n' + 1 and n' < n. Hence the updated p\ will be monic in this case. And a 
fortiori if n > 2 and s is either (i) geometric or (ii) essential and e > 1 . □ 


Example 4.15 Let n > 2 and s = sq, ..., s±- n be a geometric sequence over F with 
common ratio r = s_i/sq £ F. Proposition ^ ,1J\ yields iteratioiis ( x , so) and (x — r, s 0 )- 


Example 4.16 (Cf. Example \S.A ) Let Si = S_ t where S_ is as in Example \S.A As in 
s = a 5 , a 9 , a 4 , 0, 0, a 2 . Normalising Algorithm^ -12\ gives the following table: 


s 

A, 

a; 

e s 

T 

— 

— 

1 

1 

(1.0) 

a 5 

a 5 

1 

1 

(x, a 5 ) 

5 9 

a , cr 

a 9 

a 5 

0 

(x + a 4 , a 5 ) 

a 5 , a 9 , a 4 

a 11 

a 5 

1 

(x 2 + a 4 x + a b , a 5 x ) 

a 5 , a 9 , a 4 , 0 

a 2 

a 11 

0 

(x 2 + cP 2 a; + a 7 , a 5 x + a 11 ) 

a 5 , a 9 , a 4 , 0, 0 

a 11 

a 11 

1 

(x 3 + a 12 x 2 + a 9 x + a 4 , a 5 x 2 + a; 11 ^ + a 5 ) 

a 5 , a 9 , a 4 , 0, 0, a 2 

0 

a 11 

0 

(x 3 + a 6 a; 2 + a 3 x + a 13 , a 5 a; 2 + a 2 x + a 10 ). 


For r = s 0 ^ 0, p'{r) = p'(r,s_i) = (1,0); for — 5 < j < —2 and r = s 0 ,...,Sj, 
p!{r) = p'{r,Sj ) = /i(s 0 , • • • ,s 2 +j). Here Li = L 2 = 1, L 3 = L 4 = 2 and L 5 = L 6 = 3, 
which agrees with Theorem \ 3. 7j Note that when n — 2,4, 6 in this table and i = 1, 2, 3 in 
Example \3.2i p, = g^/lc(gi ). 
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4.4 A Worst-Case Analysis 

Next we give a worst-case analysis of Algorithm 14.121 For n > 1 and s = s 0 , • • •, s i_ n 
define a n = J2i=o L( s o,..., Sj). The following identity and inequality seem to be new. 

Proposition 4.17 If s = so ,..., si_ n then <r n = L n (n + 1 — L„) < (n + l) 2 /4, with 
equality if and only if n = 2L n — 1. 

PROOF. The equality is trivially true if s = 0 n . For the sequence s = 0" _ 1 ,s v with 
—v > 0, we have cr n = n = n(n + 1 — n) as required. Suppose inductively that n > 2 , s 
is non-trivial, the equality is true for s — s 0 ,..., s\_ n and t = s, s_ n . If A(pp, t) — 0 then 
L „ +1 = L n and by the inductive hypothesis 

<r n _|_i c n 4- L n L n (n -f 1 L n) -\- L n L n (n 4-2 L n ) L n _^i(n 4-2 

If A(/ix; t) 7 ^ 0 we apply Lemma 14.21 or 14.71 If n + 1 — L n < L n then L n+1 = L n and we 
have just seen that the result is true in this case. If n-|-l —L n > L n then L n+1 = n-|-l — L„ 
and by the inductive hypothesis, 

cr n +i = cr n + {n 4- 1 — L n ) = L n (n 4- 1 — L n ) 4- (n + 1 — L n ) = (L n 4- l)(n 4- 1 — L n ) 

Secondly, the right-hand side is (n + 1 — L n )(n + 2 — (n+1 — L n )) = {n + 1 — L„)(L„, 4 -1) 
which we have just seen is a n+ 1 . This completes the inductive proof of equality. 

For the inequality, we show that 4L n (n -f 1 — L n ) < (n + l) 2 . For integers a, b we have 
4 ab < (a + b) 2 , with equality if and only if a = b. Put a = L n and b — n + 1 — L n . Then 
a + b = n + 1, so that 4 ab < (n 4-1) 2 , with equality if and only if L n = n 4~ 1 — L n . □ 


Corollary 4.18 
linear in n 


(Of. fjfj) Let s = s 0 ,..., Si_ n be a sequence over D. Ignoring terms 
the number of multiplications in Algorithm f.12 to compute /i\ G MP(s) or a 


minimal solution p for s is at most cn 2 /4 where c is given by 


D 

outputs 

c 

domain 

hi 

3 

domain 

solution /i 

5 

field 

monic /q 

2 

field 

solution p, with monic p± 

3. 


Proof. For 1 < k < n — 1, let r = s 0 ,...,si_ fc and p\ G MP(r). Then A (yupr, s_ fc ) 
requires at most L*, 4 - 1 multiplications and G MP(r,«_(,) requires at most L;,. + 1 if 
r is pseudo-geometric and L*. 4- L y + 2 otherwise. If r is essential then by construction 
L*;/ < Lfc/ + i = • • • = Lfc so that v\ requires at most 3Lfc-|-2 multiplications. Thus computing 
a minimal polynomial for s = So,..., Si_ n requires at most Xlfc=i(3Lfc 4 - 2) < 3?i 2 /4 4- 2 n 
multiplications by Proposition 14.171 If (pi,xp 2 ) is a solution for r then \p 2 \ < L fc — 1 
and /j / 2 = 0 or \p' 2 \ < L — 1, so that we need at most L fc 4- L fc / additional multiplications 
to obtain v 2 . Ignoring linear terms, this gives at most 5n 2 /4 multiplications to obtain a 
solution for s. The remaining cases are similar. □ 
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4.5 An Identity for p and p! 


We prove an identity satisfied by p,p'. This is our analogue of Identity (00) satisfied by 
partial quotients; see Theorem 13 .1 1 First a non-zero scalar: 


Definition 4.19 We define V s G D 
initialisation. Let p\ G MP(s) and t = 


x using Algorithm \4-12\ as follows: V. s = 1 on 
s,a. If Ai = A(/ii ; f) = 0 put V t = V s ; otherwise 


f A[ V s if e s < 0 
0 Ai V s otherwise. 


If D is the field of two elements then V s = 1 for any s. Suppose that s = a,b,c with 
a, b G D x and c G D as in Example 14.31 After the first iteration, p = ( x,a ), pi = (1,0) 
and V a = a. Next p = ( ax — 6, a 2 ), p' = (1, 0) and F7 a ,b = A X V a = a 2 since e a ,b = 0. If 
Ai = ac — b 2 0 then e a ,b,c = 1 and V a ,b,c = Ai= (ac — b 2 )a 2 . 


Proposition 4.20 (Cf. |3 Theorem 7.42]) If p, pi are as in Algorithm ^. 12\ then 

p 2 hi hi h 2 V s . 


PROOF. If s is trivial, p = (1, 0), p' — (0, —1) and p-zhf — P 1 P 2 = 0 • 0 — 1 • (—1) = 1. 

Suppose inductively that P 2 P 1 — Pi p '2 — V s and t — s,a. If A x = 0, there is nothing to 

prove. Otherwise let e s <0. By construction v' = p! and 

v 2 v[ - vi u ' 2 = (A [p 2 - A ix~ e p' 2 ) u[ - (A)//, - Aix~ e p[) u ' 2 

= (A [p 2 - Aix~ e p' 2 ) p\ - (A[pi - Ai x~ e p[) p 2 
= A x (P2P1 hih^) = ^1 A 7 s = 


whereas if e s > 1 we have z/ = p and by construction 

u 2 - ui u 2 = (A [x +e p 2 - Ai p 2 ) u[ - (A [x +e pi - Ai p[) u 2 
= (A [x +e p 2 ~ Ai/4) Pi - (A[x +e p 1 - Axp’f) p 2 
= Ai (p2p[ - P 1 P 2 ) = A X V S = Vt. □ 

If A x = ac — b 2 0 in Example I4.3f ii). we have seen that 'V a ,b,c = a 2 A\ and 
P 2 P 1 — P 1 P 2 = afi’x (ax — b) — ( a 2 x 2 — abx — A x )a 2 = a 2 A x . 

We have the following immediate consequence of Proposition 14.201 

Corollary 4.21 If s is a finite sequence over F then gcd (pi,p 2 ) = gcd (pi, p[) = 

gcd(/i 2 , p'2) = I- 

The next useful consequence of Proposition 14.201 is worth stating separately. The proof 
is similar to that of Proposition 13.111 and is omitted. 

Proposition 4.22 Let f G R 2 . If rn — f 2 p. 1 — fi /o and ml = / 2 /i) — f\p 2 then 


V s /i = m'pi - rn p\. 


Example 4.23 For s — 0 n 1 ,a as in Example \2.If( ii), p 
V s = a. For f G R x , m = f 2 pi - j\p 2 = f 2 ^ - fia, ml = 


= (x n ,a), pi 
f 2 hl ~ flh 2 = 


= (1,0) and 
/ 2 and 


m'pi - mp'i = f 2 x n - (. f 2 x n - fi a) = a fi = V s fi. 
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5 Decomposition 

We now turn to the set of annihilating polynomials of a finite sequence s over D (which 
may be /S'|rz for some infinite sequence S over a held). 

We will characterise the annihilating polynomials which uses a pairing R 2 x R 2 —y R. 
This pairing was suggested by Identities (j2J) and ([3]) of the Introduction. Even though 
our conclusions for pseudo-geometric sequences turn out to be a special case of those for 
essential sequences, we have treated each case separately as their proofs differ, and little 
would be gained by combining their proofs in one place. Moreover the simpler pseudo¬ 
geometric case acts as a precursor to the remaining case. For essential sequences, the 
integer n' and the identity L n + L n > = n' + 1 are vital. In each case, we characterise the 
elements of Ann(s) using the pairing and show that if we restrict to annihilators of degree 
at most n, our decomposition is unique and we can describe the set of solutions. 

These proofs are valid once we know either a minimal polynomial or a ’minimal sys¬ 
tem’ (see Definition 15.8jl for the original finite sequence i.e. they do not depend on the 
provenance of the minimal polynomial. 

5.1 A Pairing 

Propositions 13.111 and 14.201 suggest the following definition: 

DEFINITION 5.1 For a sequence s we define a pairing ( , ) = ( , ) s : R 2 x R 2 —)• R by 

(f, 9 ) = f‘2 9\ ~ fi 9i 

where x f 2 = [fi ■ s] and similarly for g 2 . 

For s,g,t and h as in Lemma [2.61 the proof of Lemma [2.61 shows that (g, h)t 7 ^ 0 and 
\{g,h) t \ = \g\ + \h\ - n- 1 > 0. 

From Corollary 13. 101 and Proposition 13. Ill we have (—= (/, q')qi — (/, q) q[. We 
can restate Proposition 14.201 as (//, p!) = V s and Proposition 14.221 as 

VJi = {f,p')p 1 - {f,p) p[- 

5.2 Geometric Sequences. II 

Throughout this subsection, n > 1 and s — s 0 , ■ ■ ■, si_ n is a pseudo-geometric sequence 
over D. We assume that A G R 2 is a minimal solution for s, L n = • • • = Lx = |Ai| = 1 and 
A' = (1, 0), so that A 2 = (A, A') = V G D x . For example, if A is obtained via Proposition 
13.51 then A = ((x — r)/S 0 , 1) and V = 1. If A is obtained from Theorem 14.91 then X\ — x 
and V = so if s_i = 0; otherwise Ai = sgx — s_i and V = Sq. I 11 both cases we have 
V/i = (/, A')Ai — (/, A)A' 1 . 
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5.2.1 Annihilating Polynomials 

Lemma 5.2 Let f\ G Ann(s) x and x f 2 — [fi ■ s]. If (p G R, gi = fi — p>, x g 2 = [gi ■ s] 
and \ip\ < |/i| - n t/ien (%> |#i| = |/i|, 5-1 G Ann(s) x and g 2 = f 2 ; (ii) (g, f) = ipf 2 . 

Proof, (i) Firstly \ip\ < |/i| — n < |/i| — 1, so \gi\ = |/i|. Since /1 G Ann(s) we can 

write fi ■ s — F + x f 2 where v(F) < |/i| — n = |< 7 i| — n. Then 

gi ■ s = (/1 - <?) • S = F + x f 2 - <p ■ s 

and v(<p-s) = \<p\ +v = \ip\ < \fi\ —n= \gi\ —n since Li = 1 implies that v = 0. Therefore 
v (gi • s — xf 2 ) < |< 71 1 —n, gi G Ann(s) x and r / 2 = / 2 . (ii) From Part (i) and the definition 
/) = 92)1 ~ gif 2 = / 2/1 - (/1 - <f)f 2 = <pf 2 - □ 

Proposition 5.3 Let f g R 2 and x f 2 = \ ■ si. If m = If, X) and m! = (f, A') then 

m = f 2 Xi - V / 1 , m' = f 2 and \ f 2 \ = |/i| - 1. Further 

(i) f\ G Ann(s) x if and only if \m'\ — |/i| — 1 and \m\ < |/i| — n; 

(ii) if fi G Ann(s) x and \ fi\ < n then m G D. 

Proof. Since A 2 = V, Ai • s = M + x V where v(M) < 1 — n < 0. (i) For any fi G R x , 
m = {f,ty = / 2 A 1 -/ 1 V and wl = (/, A') = / 2 . We have |m'| = |/ 2 | = v+l/il — 1 = |/i|-l. 
Let fi G Ann(s) x so that fi ■ s = F + x f 2 where v(F) < \f±\ — n. Then 

fi • (M + x V) = /1 • (Ai • s) = Ai • (/1 • s) = Ai ■ (F + x f 2 ) 

xm = x (/ 2 Ai — V / 1 ) = fi-M — Xi-F and |m| + l < max{v (/ 1 • M), v(Ai-F)} < |/i| + l — n 

as claimed. Conversely, suppose that \m'\ = |/i| — 1 and \m\ < |/i| — n. We claim that 

m'Ai G Ann(s): 

(m'Ai) ■ s — ml ■ (M + Vi) = m' ■ M + x V m' 

and v(m / • M) < |m'| + 1 — n = \mXi\ — n. We have m = / 2 Ai — V/i = m!X\ — V/i, so 
V/i = m'Ai — m. Consider — V fi: \gi\ = |/i| and |m| < |/i| — n = \m'\ + 1 — n = 
|m'Ai| — n since |m| < |/i| — n = |^i| — n. So g\ G Ann(s) by Lemma 15.21 and hence 
fi G Ann(s) x . (ii) This is immediate. □ 


Corollary 5.4 (i) V Ann(s) x C {ip'X 1 — p : p>' ^ 0, \<p\ < \ip'\ + 1 — n} c Ann(s); 
(ii) VMP(s) C {ip'X 1 -ip : ip'E. D x , | V |<l-n}CMP(s); 

(in) if n > 2 then VMP(s) C {ip'X 1 : ip' G D x } C MP(s). 

Proof, (i) If f\ G Ann(s) x then V f\ — m'X 1 — m where \m'\ = |/ 2 | = |/i| — 1 > 0 and 
\m\ < | f\\ — n — \m'\ + 1 — n by Proposition 15.31 If p' 7 ^ 0 and |(^| < \ip'\ + 1 — n then 
ip'X 1 — ip E Ann(s) x by Lemma [5.21 (ii), (iii) These are immediate. □ 
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5.2.2 Solutions 

We apply the results of the previous subsection to finding solutions for a pseudo-geometric 
sequence; this is a precursor to the discussion of solutions for essential sequences in Sub¬ 
section 15.3.21 

Lemma 5.5 Let p,p' G R. If g i = p'Xi — p, 0 < \p'\ < n — 1, p G D and xg 2 — [g i ■ s] 

then g 2 = p'X 2 — V p'. V p = (g, A) and V p' — (g, A'). 

Proof. We have Ai • s = M + x X 2 where v(M) < 1 — n and 

gi ■ s = p' ■ (Ai • s) — Lp ■ s = <p' ■ (M + x A 2 ) — <p ■ s = (p' ■ M — tp ■ s + x ip'X 2 . 

Further, ■ M) < |<p'| + 1 — n < 0 by hypothesis and ■ s) = v = 0 since i^eD. 
Thus v{ip' ■ M — ip ■ s) < 0 and g 2 = ip'X 2 = V ip'. We have 

(g, A) = .92 Ai - g 2 X 2 = V ip'Xi - ( p' Ai - ip )V = V <p 

and (g, A') = g 2 X[ - g^ = g 2 = V ip'. □ 


Corollary 5.6 Let fi e Ann(s), |/i| < n and m! = (/, A'), m = (/, A). Then 

(%) (uniqueness) if V f\ = ip'Xi — p where 0 < \ip'\ < n — 1 and p e D fhen p — m 
and p' = m'; 

(ii) V / 2 = m'A 2 — mA' 2 . 

Proof, (i) Applying Lemma l5Al to (?i = V /i = <p'Ai — p gives V p = (g, A) = V(/, A) = 
Vrn. Hence p = m = f 2 and similarly p' — ml — f 2 X\ — V f±. (ii) We have V / 2 = 
m'A 2 — mAf, since m' = / 2 , A 2 = V and A 2 = 0. □ 

Thus if s is a geometric sequence over F, p, V are obtained from Algorithm 14.121 and 
// = (1,0) then the minimal solutions for s are 

{(tV-vW) : e F x , \p\ < 1 -n}. 

Also if Sq 7 ^ 0 and n G [ 1 , 712 ) we may take A = (ai, 1), A' = (1,0) and V = 1. For 
example, from Corollary 15.41 

Corollary 5.7 Let S be an infinite geometric sequence over F, r = S-i/So, n G [l,n 2 ) 
and s = S\n. The minimal solutions for s are 


{(p'(x-r)-p,p'): p' G F x , \p\ < 1 -n}. 










5.3 Essential Sequences. II 

When s is essential, more information is available for decomposition. Informally, we have 
a pair of linked triples, their second components and V 6 D x , all related by the pairing 
of Definition 15.11 

Definition 5.8 Let n > 2 and s = s 0 , be an essential sequence over D. A 

minimal system for s is a 5-tuple (A, n', A', (, ) s , V) consisting of 

(i) a minimal solution A G R 2 for s and Ai fL Ann(s, s_ n ); 

(ii) n' = maxi< j<n {j : Lj < L n } and s' = s 0 ,..., si- n >; 

(Hi) a minimal solution A' G R 2 for s', X[ fL Ann (s', s_ n /) and L n + L n > — n! + 1; 

(iv) the pairing {,) s :R 2 xR 2 gRo/ Definition ^. 1[ 

(i;jV=(A,A') s GDA 

From Theorem 13.II and Corollary 13. 101 (q, Hi — 1, q', (, )s| n , (—l)*^ 1 ) is a minimal sys¬ 
tem for S\n if n G [ni,rii + 1 ) and S\n is essential. We have seen that if (n, g i,A) and 
(ra',/A, A') are linked triples for s then (//, n',g', (, } s , V s ) is a minimal system for s by 
Theorem 14.91 Definition 14.191 and Proposition 14.201 

N.B. Throughout this subsection, n > 2, s is a sequence over D and (A, n!, A', (, ) s , V) 
is a minimal system for s — s 0 ,..., Si_ n . We put ( , ) = ( , ) s , L = L n and L' = L n /. 

As we have already seen in Propositions 13. 1 fl and l4~20l ii) . for any / G R 2 we have 

V /i = (/, A / )A 1 — (/, A)A , 1 . 

5.3.1 Annihilating Polynomials 

Lemma 5.9 If fi G Ann(s ) x 7 G R, g\ — fi — and \<p\ < |/i| + L — n — 1 then 
(i) \9i\ = \fi\, 9i e Ann(s) x and g 2 = f 2 ~ ^A' 2 ; 

(a) (g,f) = <p{f,y)- 

Proof, (i) Firstly |gi| = |/i| since n' < n implies that 

I^A'il = |v?| + L' < |/i| + L — n — 1 + L' = |/i| -n + n < |/i| - 1. 

Since fi G Ann(s) we can write f\ ■ s — F + x f 2 where v(F) < \ f \ \ — n — \gi\ — n and 
A'[ - s' = M' + x X' 2 where v(M') < L' — n' = 1 — L. Put N' = X\ ■ (s — s'). Then 
v(N') < L' — n' = 1 — L and so A • s = N' + M' + x A ' 2 where v(N' + M') < 1 — L. Thus 

gi ■ s — (f 1 -(p\[)-s = F + xf 2 -ip- (N' + M' + x A 2 ) = F -ip - (A 7 + M') + x (/ 2 - 2 ). 

Now v((y9•(A / -|-M , )) < |<£>| + 1—L < |/i|+L— n —1 + 1—L = \fi\—n— \gi\— n, g± G Ann(s) x 
and g 2 = f 2 - p\ ' 2 . (ii) From Part (i) ( g , /) = (/ 2 - 97 A 2 )/i - (A - v?A \)f 2 = p{f , A'). □ 
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Theorem 5.10 Let f G R 2 . If m = (f, A) and ml = (/, A') i/ien 

(i) /i G Ann(s) x if and only if |m'| = |/i| — L and \m\ < |/i| + L — n — 1; 

(ii) if fi £ Ann(s) x t/ien |m| + L' < |/i| — 1. If in addition \fi\ < n then |m'| < n — L 
and |m| < L — 1. 

Proof. First write Ai • s — M + x A 2 where v(M) < L — n < 0. (i) Let f\ G Ann(s) x , 
so that /1 • s — F + x f 2 where v(F) < |/i| — n. Then 

fi ■ (M + x A 2 ) = fi ■ (A 1 • s) — A 1 • (/1 • s) = A 1 • (F + a; ff) 

xm — x (/ 2 A 1 — /iA 2 ) — fi ■ M — Xi ■ F and |m| + 1 < max{|/i| + v(M), L + v(F)} < 
|/i | + L — n. Hence 

IH + L' < |/i| + L - n - 1 + L' = |/i| +n' - n < \fi\ - 1 

as nl < n. We have V fi = m'A 1 — mA'j, so |m'| = |/i| — L. Conversely, if |m'| = |/i| — L 
then m'A 1 G Ann(s), for 

(m'A i) • s — m! ■ (M + x A 2 ) = m' • M + xm'X 2 

and v(m' • M) < |m'| + L — n = |m'A i| — n. We claim that V /1 = m'A 1 — mA ' x G Ann(s). 
From Lemma [5.91 it suffices to check that |m| < |m'Ai| + L — n — 1. But |m'| = |/i| — L, 
so |m| < |/i| + L — n — 1 suffices, and this is true by hypothesis. We conclude that 
V/i G Ann(s) and hence so is f\. (ii) The first sentence was proved in Part (i); we also 
have |m'| = |/i| — L < n — L and |m| < |/i| + L — n — 1 < L — 1. □ 

Recall that for any sequence s — so, • • •, s i-n over D, e s — n + 1 — 2L n G Z. 

Corollary 5.11 (Cf. [7J) 

(i) V Ann(s) x C 1 — (^A^ : 7 ^ 0, |<p|<|<p'|—e s }CAnn(s) x ; 

(njVMP(s)C{^Ai-^Ai : /eD x , |^| < -e s } C MP(s); 

(Hi) if 2L < n then VMP(s) C : (p f G D x } C MP(s). 

Moreover if V is a unit of D (/or example if D is a field) the inclusions are equalities. 

Proof, (i) If j\ G Ann(s) x then V/i = m'A 1 — rnX\ where |m'| = |/i| — L and 
|m| < |/i| + L — n — 1 = |m'| — e s by Theorem 15.101 If 99 ' 7 ^ 0 and |<^| < \<p'\ — e s 

then (p'X 1 — (pA ( G Arm(s) x by Lemma 15.91 (ii) If /1 G Ann(.s) and |/i| = L then by 

Theorem 15.101 |m'| = 0. If /1 = (p'Ai — (pA), <p' G D x and |<p| < — e s then fi G Ann(s) 
and |/i| = |Ai| = L by Lemma 15.91 Part (iii) is an immediate consequence of (ii). If 
fi G Ann(s) x and V is a unit of D then /i/V G Ann(s) x and hence f± G VAnn(s) x . □ 

Thus if e s < 0 we have Ai — A' x G MP(s), as is well-known for sequences over a field. From 
Theorem 13.71 and Corollary 15.111 we have 
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Corollary 5.12 (Cf. JE, Theorem 1]) Let S be an infinite sequence over F, n G 
[rii,n i+ 1 ) and q = q^ l \ q' = gh-R If s = S\n is essential then 

Ann(s) x = {p'qi-pq[ : p' fi 0, \p\ < |y/|-e a },MP(s) = {( p'qi~pq{ : <p' G F x , |^| < -e s }. 

5.3.2 Solutions 

Next we look at solutions i.e. pairs (/i,/ 2 ) with /i G Ann(s) x and x f 2 — [fi ■ s]. 

Lemma 5.13 Let p, ip' G R. If g\ = </?'A i — <^A ' l7 0 < \ip'\ < n — L and \ip\ < L — 1 then 

(i) g 2 — ip'X 2 — 99 A 2 ; (ii) V p = (g, A) and V ip' = (g, A'). 

Proof, (i) We have Ai-s = M + xA 2 and A ' x ■ s' = M' + a; A 2 where v(M) < L — n 

and v(ilf') < L' — n' — 1 — L. Write s = (s — s') + s' so that v(s — s') < —n' and put 

N' — A [ ■ (s — s'). This gives 

fi'i-s = (<£>'Ai —^Ai) • s = ip' ■ (M + x\ 2 ) - ifX[ ■ ((s - s') +s') 

= ip' ■ (M + x \ 2 ) — ip ■ N' — ip ■ (. M' + x A 2) 

= ip' ■ M — ip ■ N' — ip ■ M' + x {ip'X 2 — <p\ 2). 

Further, v(<^' • M) < |(^'| + L — n < 0 by hypothesis and similarly v(<£> • M') < 0. Now 
v(<^ • N') = \ip\ + lA'J + v(s - s') < M + L' — n' — \<p\ + 1 — L < 0 
as \ip\ < L — 1. Thus v(<// • M — ip ■ N' — ip ■ M') < 0 and g 2 = ip'X 2 — tp\ ' 2 . (ii) We have 

{d: ^) — ^2 A 1 — (?i A 2 = (ip' A 2 — ip A 2 ) A 1 — (ip' A 1 — ip A',) A 2 
= ip (A 2 A) — A 1 A 2 ) = V</9. 

Similarly (g, X') = g 2 A 1 — pi A 2 = 2 A^ — A 1 A ' 2 ) = V</?'. □ 

Corollary 5.14 Let f be a solution for s, \fi\ < n and ml — (/, A') ; m = (/, A). 

(%) (uniqueness) If V /1 = <^'A 1 — (^A^ where 0 < |<^'| = |/i| — L and |</?| < L — 1 then 
ip = m and ip' = m' ; 

(ii) V / 2 = m'A 2 — mA ' 2 ; 

(degree bound) if X' 2 fi 0 then \m\ + |A 2 | < \f 2 \ — 1 and \f 2 \ = \m'\ + | A 2 | - 

Proof, (i) Applying Lemma 15.131 to g\ — V f\ — p'X 1 — pX [ gives V m = V(/, A ) = 

(V f, X) — (g, X) — V <p. Therefore p — rn and p' = ml. 

(ii) We know from Theorem 15.101 that V fi — m' X 1 — rn X \ where |m'| < n — L and 
\m\ < L — 1. Hence Lemma [5.131 implies that V / 2 = m!X 2 — mXfi 

(iii) We have A' G Ann(s') and if s' is trivial then A' = (c, 0) for some c G D x . As 
A 2 fi 0, s' is non-trivial so 1 — n' < v(s') < 0 and therefore v(s') = v. From Theorem 
15.101 |m| + L' < |/i | - 1, so 

M + IA 2 I = \m\ + (v(s') + (A'J — 1) = \m\ + v + L' — 1 < |/i| — 1 + v — 1 = |/ 2 | — 1. 

From Part (ii) we have V / 2 = m'A 2 — mA' 2 , so |/ 2 | = |m'| + |A 2 1- Cl 

The final result of this section on solutions is a simple consequence of Corollary 15.141 
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Corollary 5.15 (Cf. w If n > 2, s — sq, ..., si- n is an essential sequence and £ 
denotes the solutions {(/i, f 2 ) '■ fi G Ann(s), 0 < |/i| < n} then 

V £ C {<£>' A — <p A', 0 < \tp'\ < n — L, \tp\ < \tp'\ — e s } C £ 

and ifV is a unit of D (for example if D is a field) the inclusions are equalities. 

We leave the corresponding result for minimal solutions to the interested reader. 

6 Some Applications of Decomposition 

We give some applications of the results from the previous sections. As usual, n > 1, 
s = s o, • • •, Si_ n is non-trivial and p, p' are obtained using Algorithm 14.121 or. if D is a 
field F, using the Normalised Algorithm 14.121 We put L = L(s). 

6.1 Sequences over a Field 

We prove several gcd-related results, relate partial quotients to p i, p\ and count the 
number of solutions when |F| < oo. Firstly a partial converse to Proposition 12.8( h). 

Corollary 6.1 7/2L < n, f\ £ Ann(s) x and \fi\ <n — L then 

(i) V s / = m'p; 

(ii) if in addition gcd(/i,/ 2 ) = 1 then m! £ F x i.e. f\ £ MP(s). 

Proof, (i) Since L < n — L, such an f\ can exist. Proposition 15. 31 or Theorem 15.101 imply 
that V s /i = m' pi — mPi where \m'\ = |/i| — L > 0 and \m\ < |/i| + L — n — 1 < —1, 
so V s /i = m'p. Since |/i| < n, we also have V s / 2 = ml p 2 — m p' 2 from Corollary 15.41 
or Corollary 15.141 i.e. V s / = m'p. (ii) By Corollary 14.211 gcd(pi,p 2 ) = 1, so V s = 
gcd(V. s /i, V .5 f 2 ) = gcd (m'pi,m'p 2 ) = ml. Thus m! £ F x and |/i| = |pi| = L. □ 

The example after Proposition 12.81 shows that the condition |/i| < n — L is necessary. 
Secondly, one may show directly that if /i, pi £ Ann(s) x and |/i| + |gi| < n then (/, g) = 0; 
see pm Corollary 3.25]. This gives another proof of Corollary 16.11 

Corollary 6.2 (Cf. p. 439-444])- LC S be a linear recurring sequence over F, 
Ids = 9i F[x] where g\ is monic, n > 2\gi\ and s = Ain. Then 

(i) g is a minimal solution for s and any minimal solution of s is eg for some c £ F x ; 

(ii) S_ — [<71 • s]/< 7 i and S is determined by S 0 , ..., A 2 | 5l |_i; 

(in) if i is the first index such that b t+ i = 0 in obtaining the partial quotients of S_, 
then g = q^/lc(q[ l ' > ). 

Proof, (i) Firstly, g is a solution for s = S\n. As is a minimal polynomial of S, 
gcd(< 7 i, g 2 ) = 1 and xg 2 = [rp-A] = [r/, -s] by Lemma E31 since n > Ig^. As |^i| < n-\gi\, g ] 
is a minimal solution for s by Corollary 16.11 Further e s > 0 so any minimal solution for s is 
eg where c £ F x by Corollary 15. 41 or I5TI1 (ii) We have A = [grS\/gi = [<7i-s]/<7i = xg 2 /g\. 
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Since g\ is uniquely determined by so,..., s 2 | si |_i, so are g 2 and S. (iii) If s is geometric, 
this is Example 14.151 Suppose that s is essential and put q = q^\ q' = q^~ l \ We have 
n i+ 1 = oo since S_ G F(x) and n* = \q[\ + |gi| < 2\qi\ = 2L < n < n t+ 1 . From Theorem 
13.71 q is a minimal solution for s and g is the unique monic solution of s by Corollary 15.111 
since e s > 0 and L < 2L < n. □ 

Instances of Part (iii) of Corollary 16.21 were given in Example 14.161 

Corollary 6.3 Suppose that f is a solution for s such that \ f\ \ < n and let m = (/, p), 
ml = (f,p r ). Then gcd (m,m/) = gcd(/i,/ 2 ). 

PROOF. By definition, m = / 2 /o — f\ /i 2 and ml = f-ipf — f\pl 2 so that if d |/i, / 2 then 
We also know that V/i = mlp\ — mp' x by Proposition 14.201 Corollary 15.141 
implies that V / 2 = mlp^ — mp' 2 since |/i| < n. Hence if d \m, ml then d |/i, / 2 . □ 

For the next result, ¥ q is a finite field with q < oo elements. If d > 0, the number 
of polynomials with coefficients in ¥ q of degree d is — (q — 1 )q d and the number of 
polynomials of degree at most d is 1 + Ylt =o Results of Section [5] now easily give the 

number of solutions for s with denominator of degree d when L < d < n: 

Corollary 6.4 For L < d < n, the number of solutions for s with denominator of 
degree d is N d _ L (l + ■ 

PROOF. Let E d = {(/i,/ 2 ) : /i e Ann(s) x : |/i| = d}. From Corollary 15.61 or Corollary 
15.111 we have / G E d if and only if /i = tp'pi — ipp\ where (i) \ip'\ = d — L and (ii) tp = 0 
or \<p\ < d — L — e s , which yields the stated result. □ 


6.2 Non-Vanishing Annihilating Polynomials 

We consider the following problem: let a G D be arbitrary and suppose that pi(a) = 0. 
Find a solution £ = (£i,£ 2 ) such that £i(a) ^ 0 and £i has least degree among solutions 
with first component not vanishing at a. We begin with a pseudo-geometric example. 


Then e s = n — 1 > 0 and 
Thus all minimal polynomials 


Example 6.5 Let n > 2 and s = sq, . ..,si_ n = l,0 n_1 . 

V s = 1, so MP(s) = {</?'x : <p' G D x } by Corollary \5.f 
of s vanish at 0. However g = ( x n + l,x n_1 ) is a solution for s and pi(0) 7 ^ 0. We will 
shortly see that min{|/i| : f is a solution for s, /i( 0 ) 7 ^ 0 } = n, so that |pi| attains this 
minimum. 


We can assume that s is non-trivial and L > 1, for otherwise (Ji G D x vanishes 
nowhere. Put Ann(s)^ = {/1 G Ann(s) : fi(a) 7 ^ 0}. Any polynomial of degree n which 
does not vanish at a annihilates s, so that 

= mind/! | : fi G Ann(s)^ a -*} 
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is well-defined and L < id 0 ') < n. We put MP(s)^ = {/i G Ann(s)) a ) : |/i| = L( a )}. 
Example 16.51 shows that L' 0) — L can be arbitrarily large. As usual, // is obtained as in 
Algorithm 14.121 

Corollary 6.6 If n\(a) = 0 then /A(a) ^ 0. 

PROOF. Proposition 14.201 yields fi 2 fA — // 1 /i 2 G F x , so fi\ (a) ^ 0 (and /i 2 (a) A 0). □ 

Using Proposition 15.31 Theorem 15. 101 and Corollary 16.61 we can now solve the problem 
posed at the head of this subsection. 

Theorem 6.7 (Cf. [lf\ Proof of Theorem 3.7] ) Let n > 1, so,..., si_ n he a sequence 
over D, e = n + 1 - 2L n and M = max{e, 0}. If //i(a) = 0 then ld“) = L + M. In fact 
£i = x M /ai — ji\ G MP(s)W and £ 2 = x M [i 2 — n' 2 . 

Proof. If n = 1 and Hi( a ) = 0 then a = 0 and e = 0; £i = Hi — Hi = x — 1 G MP(s) 
satisfies £i(a) A 0 and £2 = So = H 2 — /x 2 , so £ is the required solution and ld a) = L + M. 

If n > 2 and e < 0 then £1 — hi — lA £ MP(s) so L = L( a ) and £ 2 = H 2 ~ lA by Lemma 
15.131 Also £i(a) A 0 by Corollary 16.61 and so £ is the required solution. 

Now let n > 2 and e > 0. We have L + e = n + l — L < n since L > 1. We first show 
that L ( °) > L + e. Let f\ G Ann(s)^, / = (/i,/ 2 ), m = (/, n) and m' = (/,//). From 
Proposition 15.31 or Theorem 15. 101 we have V s f\ = mln 1 — where \m'\ — |/i| — L > 0 

and \m\ < |/i|+L— n— 1. If |/i|+L— n —1 < 0 then m = 0, V s /i(a) = m'(a ) Hi (a) = 0 and 
fi A Ann(s)^ for a contradiction. Hence |/i| > n+1 —L = L+e and L( a ) > L+e = L +M. 

To see that L^ < L + M, let £1 = x M Hi — h\ which has degree L + M. We have 
£i(a) = —Hi(a) 7 ^ 0 by Corollary 16.61 We claim that £1 G MP(s)(“). We have |1| = 
0 = |£i| + L — n — 1, so by Lemma 15.21 or Lemma 15.91 we have £1 G Ann(s , ) l - a ) and 
L (a) < |£i| = L + M. 

Finally, we verify that £ 2 = x M H 2 ~ H 2 - Put T' — XjM and T = 1- If s is pseudo¬ 
geometric, Lemma 15.51 applies since |<p| = 0 < L — 1 and 0 < M < n — 1. Hence 
£ 2 = x M H 2 — H 2 = x M H 2 - Suppose that s is essential. We have 0 < M < n + 1 —2L < n — L 
since L > 1, so that Lemma 15.131 applies and £ 2 = x m /i 2 — /i 2 in this case too. □ 

Theorem 16.71 yields the following simple extension of Algorithm 14.121 

Algorithm 6.8 (Cf. [TQ Algorithm 3.2]) 

Input: n > 1, a G D and sequence s — sq, ..., s 1 _ n over D. 

Output: Solution £ for s such that £1 G MP(s) 6 f). 

|~ Algorithm 14.121 (input : n, s; output : HiH')‘i 

if Hi{°) A 0 then £ «— /1 else £ G- X max ^ n+1 ~ 2 ^’ 0 ^ h ~ h'i 

return £.J 
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Tabic 2: Algorithm 14.121 for 0,1,1, 0, 0,1, 0,1 over F 2 


s 

A x 

e-s 

T 

s’ 


— 

— 

( 1 . 0 ) 

( 0 , 1 ) 

0 

0 

1 

( 1 . 0 ) 

( 0 . 1 ) 

0,1 

1 

2 

(^, 1 ) 

( 1 . 0 ) 

0 , 1,1 

1 

-1 

( X 2 + X, 1 ) 

( 1 . 0 ) 

0 , 1 , 1,0 

1 

0 

(x 2 + X + 1 , 1 ) 

( 1 . 0 ) 

0 , 1 , 1 , 0,0 

1 

1 

(x 3 + x 2 + X + 1 , x) 

(x 2 + X + 1 , 1 ) 

0 , 1 , 1 , 0 , 0,1 

0 

0 

(x 3 + x 2 + X + 1 , x) 

(x 2 + X + 1 , 1 ) 

0 , 1 , 1 , 0 , 0 , 1,0 

1 

1 

(x 4 + X 6 + 1 , x 2 + 1 ) 

(, x 3 + x 2 + X + 1 , x) 

0 , 1 , 1 , 0 , 0 , 1 , 0,1 

1 

0 

(x 4 + x 2 + X, x 2 + X + 1 ) 

(x 3 + x 2 + X + 1 , x). 


Example 6.9 Table 2 gives the iterations of Algorithm ] f (implemented using J \3jj) for 
the sequence s over¥ 2 Table If; we have omitted as it is the constant 1. We see 

that s is essential, e = e s = 1 and fi i(0) = 0. Theorem \6. 7| implies that L (0) = L + e = 5 
and Algorithm \6.A gives £ = x e p, + p! = (A 5 + x + 1, x 3 + x 2 ). 


Remarks 6.10 (i) Algorithm \6.A is simpler than {771 Algorithm 3.2], e.g. it does not 
include tests on ji\ . It also computes /i 2 . Corollary 1 6.1 2V n) below, a version of which was 
stated without proof in m was used to justify Algorithm 3.2, loc. cit. In Algorithm 3.2, 
loc. cit. the polynomial ji\ is initialised to 1 as in fffj rather than 0; see Remark \f.l3\ 

(ii) The original motivation of let a = 0 and be the reciprocal of £ i. Since 
£i(0) 7 ^ 0, |£{| = |^i | = S say, (£{•s*)* = where j = 1 —n+5—i, and 5+1—n < i < 0 

if and only if 5 + 1 — n < j < 0. Hence if 1 < 5 < n, and the first 5 terms so ■ .., si-,5 
uniquely determine the last n — 5 terms ss ,..., si_ n if and only if 1 < n — 5 < n, and 
the last n — 5 terms s i_ n ..., ss- n uniquely determine the first 5 terms s,5_ n+ i,..., sq. 


We can also construct an element of MP(s)^ a - ) by extending s by one term. 


Corollary 6.11 (Cf. flff ) Let s = s 0 ,... //1 G MP(s) and pi (a) = 0. Suppose 

that e — e s > 1. Put t = s, s_ n where s_ n is chosen so that A(/ii; t) ^ 0. If v i is obtained 
as in Theorem A 5| then iq G MP(s)(“). 


Proof. Let A x = A(/r x ;£). Since e > 1, z/ x — A\ x e p { — A x // x G MP(£) from Theorem 
14.91 and \ui\ — n + 1 — L — by Theorem 16.71 Further, z/ x (a) = —A x /F x (a) 0 by 
Corollary 16.61 and since Ann(f) C Ann(s), z/ x G MP(s)^. □ 


For the Example of Table 2, A 9 = 1 requires Sg = 0 and we obtain u = f as before. 

Corollary 6.12 Let s — so,... , s x _ n be a sequence over F, /i, p! be as usual and m = 
(/, p) ,ml = (/, pi). If a G F, p±(a) = 0 and M = max{e s ,0} then 
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(i) Ann(s) (o) = {^Vi - ThA W\ A 0, \y\ < \ip'\ - e s , <p(a) A 0}; 

(ii) MP(s)^ = {fi G Ann(s)W : \m'\ = M, \m\ < |/i| + L — n — 1 ,m{a) A 0}; 

= {ip'fx i - (pfx[ : \(p'\ = M, \ip\ < M - e s , <p(a) A 0}. 

Proof, (i) This is a restatement of Corollary 15.41 or Corollary 15. Ill (ii) We have V/i = 
m!Hi — rriHi and j\ G Ann(s) if and only if \m'\ = \fi\ — L and \m\ < |/i| + L — n — 1. 
Also |/i| = L + M from Theorem 16.71 since fi G MP(s)^. Similarly if g\ = (p'^i ~~ '-PlA 
then gi G Ann(s) if and only if |y/| = \gi\ — L, \<h\ < \ip'\ — e s by Lemma [5.21 or Lemma 
15.91 and |gi| = L + M by Theorem 16.71 □ 

For Example 16.91 e s = 1 and so MP(s 0 ,..., s_ 7 )^ 0) = {x Hi + h'h ( x + l)hi + lA } • 


References 

[1] Berlekamp, E. R. Algebraic Coding Theory. Series in Systems Science. McGraw Hill, 
New York-Toronto, 1968. 

[2] Cheng, U. On the Continued Fraction and Berlekamp’s Algorithm. IEEE Transac¬ 
tions on Information Theory , 30:541-544, 1984. 

[3] Cocoa Team. A System for Doing Computations in Commutative Algebra. Available 
at http://cocoa.dima.unige.it, Version 5.0. 

[4] Gustavson, F. G. Analysis of the Berlekamp-Massey linear feedback shift-register 
synthesis algorithm. IBM J. Res. Dev., 20:204-212, 1976. 

[5] Lidl, R. and Niederreiter, H. Finite Fields. Encyclopedia of Mathematics and its 
Applications. Addison-Wesley, Reading, 20, 1983. 

[6] MacWilliams, F. J. and Sloane, N. J. A. The Theory of Error-Correcting Codes. 
North Holland, Amsterdam, 1977. 

[7] Massey, J. L. Shift-Register Synthesis and BCH Decoding. IEEE Trans. Inform. 
Theory , 15:122-127, 1969. 

[8] Mills, W. H. Continued Fractions and Linear Recurrences. Mathematics of Compu¬ 
tation , 29:173-180, 1975. 

[9] Niederreiter, H. Sequences with Almost Perfect Linear Complexity Profile. Advances 
in Cryptology — Eurocrypt ’87 (D. Chaum, W.L. Price, Eds.). Lecture Notes in 
Computer Science, 304:37-51, 1987. 

[10] Northcott, D. G. Injective Envelopes and Inverse Polynomials. J. London Math. 
Soc., 8:290-296, 1974. 


36 











[11] Norton, G. H. On the Minimal Realizations of a Finite Sequence. J. Symbolic 
Computation , 20:93-115, 1995. 

[12] Norton, G. H. On Shortest Linear Recurrences. J. Symbolic Computation , 27:323- 
347, 1999. 

[13] Norton, G. H. The Berlekamp-Massey Algorithm via Minimal Polynomials. 
math.ArXiv: 1001.1597, pages 1-22, 2010. 

[14] Salagean, A. An Algorithm for Computing Minimal Bidirectional Linear Recurrence 
Relations. IEEE Trans. Eifo. Theory, 55:4695-4700, 2009. 

[15] Welch, L. R. and Scholtz, R. A. Continued Fractions and Berlekamp’s Algorithm. 
IEEE Trans, on Information Theory, 46:19-27, 1979. 


37 



